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PREFACE 


The  two  studies  reported  herein  were  performed  at  Colorado  State  University,  Fort  Collins, 
CO,  and  the  University  of  Nottingham,  Nottingham,  England,  under  contract  to  the  US  Army 
Engineer  Waterways  Experiment  Station  (WES)  during  the  period  October  1 989  to  June  1 992. 
This  investigation  was  sponsored  by  the  Headquarters,  US  Army  Corps  of  Engineers 
(HQUSACE),  under  the  Flood  Control  Structures  Research  Program  as  part  of  Civil  Works 
Investigation  Work  Unit  No.  32544,  "Riprap  Toe  and  End  Section  Design,”  under  HQUSACE 
Program  Monitor,  Mr.  Tom  Munsey. 

This  investigation  was  accomplished  under  the  direction  of  Messrs.  F.  A.  Herrmann,  Jr., 
Director  of  the  Hydraulics  Laboratory  (HL),  WES;  R.  A.  Sager,  Assistant  Director  of  HL;  and  G. 
A.  Pickering,  Chief  of  the  Hydraulic  Structures  Division,  HL.  The  Contracting  Officer's 
Representative  was  Dr.  S.  T.  Maynord,  who  was  under  the  direct  supervision  of  Mr.  N.  R.  Oswalt, 
Chief  of  the  Spillways  and  Channels  Branch,  Hydraulic  Structures  Division,  HL. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was  Dr.  Robert  W.  Whalin. 
Commander  was  COL  Leonard  G.  Hassell,  EN. 
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SUMMARY 

Bank  erosion  is  a  serious  problem  to  river  engineers  concerned  with 
channel  stabilization  and  navigation.  Severe  erosion  often  occurs  at  the 
outer  bank  in  channel  bends,  where  flow  velocities  adjacent  to  the  bank 
are  elevated  due  to  the  effects  of  curvature  on  channel  flow.  Eroding 
banks  may  be  stabilized  and  protected  from  erosion  using  riprap.  When 
selecting  the  appropriate  size  of  stone  to  be  used  to  protect  a  bank  in  a 
given  bend,  it  is  necessary  to  be  able  to  predict  the  intensity  of  flow 
attack  on  the  bank.  This  may  be  represented  by  either  the  near  bank 
velocity  or  tho  boundary  shear  stress  on  the  bank.  This  report  deals  with 
the  development  of  improved  methods  to  predict  outer  bank  velocities  and 
shear  stresses.  Two  approaches  are  examined.  The  first  uses  a  statistical 
treatment  of  observed  data  from  natural  and  artificial  channels  to 
formulate  predictive  equations  for  the  ratio  of  depth  avaraged  longstream 
velocity  over  the  toe  of  the  outer  bank  and  for  the  shear  stress  in  that 
location.  The  second  tests  two  analytical  models  of  bend  flow  to  gauge 
their  accuracy  and  set  limits  to  their  applicability  in  predicting  outer 
bank  velocity. 

The  results  show  that  several  factors  appear  to  influence  outer  bank 
velocity  at  a  natural  bend.  Multivariate  equations  involving  radius  of 
curvature  to  width  ratio,  relative  bend  length,  width  to  depth  ratio, 
relative  depth  and  bank  angle  are  proposed  to  predict  the  ratio  of  outer 
bank  toe  velocity  to  average  velocity.  Simplified  equations  using  only  the 
radius  of  curvature  to  width  ratio  are  also  proposed.  The  configuration  of 
the  channel  upstream  of  the  bend  is  shown  to  be  important,  and  separate 
appproaches  are  formulated  for  bends  downstream  of  straight  and 
meandering  reaches.  For  artificial  channels  Rc/w  dominates  the  analysis, 
but  it  is  also  shown  that  the  mobility  of  the  bed  strongly  influences  the 
outer  bank  velocity  and  shear  stress. 

Model  tests  reveal  that  the  model  developed  by  Bridge  (1982) 
consistently  predicts  the  observed  outer  bank  toe  velocity  to  within  +/- 
15%.  Errors  grow  alarmingly  for  bends  with  Rc/w  values  less  than  2  and 
the  model  crashes  for  bends  with  Rc/w  <  1.  Odgaard's  (1989)  model  tended 
to  under  predict  outer  bank  velocity  by  between  5  and  40%.  This  was  the 
case  because  the  model  did  not  predict  outer  bank  scouring  in  bends  with 
bed  material  coarser  than  medium  sand.  However,  its  application  was 
limited  because  it  predicted  negative  depths  at  the  inner  bank  and  crashed 
for  long  bends.  In  contrast  to  Bridge's  model,  Odgaard's  model  remained 
stable  at  very  low  Rc/w  bends,  errors  remaining  in  the  5  to  40%  range. 

It  is  recommended  that  the  results  of  this  study  be  further  tested  and 
verified.  However,  on  the  basis  of  the  results  to  date,  the  model  developed 
by  Bridge  is  recommended  for  use  in  bends  with  Rc/w  valus  greater  than  2. 
For  very  tight  bends,  Odgaard's  model  shows  strong  potential,  but  it  must 
be  modified  to  allow  greater  mobility  and  scour  of  coarse  bed  materials. 
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MAIN  TEXT 


Introduction 

Serious  bank  erosion  often  occurs  at  the  outer  bank  in  meander  bends. 
This  erosion  is  driven  by  the  natural  tendency  of  river  meanders  to 
increase  in  amplitude  and  migrate  downstream  through  time.  The  severity 
of  flow  attack  on  the  bank  is  known  to  be  controlled  by  the  hydraulics  of 
flow  adjacent  to  the  bank  and  especially  the  propensity  for  scour  in  the 
area  of  the  bank  toe.  Conversely,  the  mechanics  of  failure  and  the 
sequence  of  events  involved  in  the  erosion,  collapse  and  basal  clean-out 
phases  of  bank  retreat  are  closely  related  to  the  engineering  properties  of 
the  bank  materials  and  the  bank  stratigraphy.  But  the  overall  rate  of 
retreat  of  the  bank  is  known  to  be  determined  by  the  capacity  of  the 
near-bank  flow  to  entrain  and  remove  slumped  bank  materials,  while 
continuing  to  erode  the  bank  and  ♦  ,gger  further  failures  (Thorne,  1982; 
Lapointe  and  Carson,  1986). 

The  importance  of  bank  attack  and  toe  scour  by  the  flow  have  long 
been  recognised,  and  their  intensity  has  been  found  to  be  a  function  of  the 
ooundary  shear  stress  acting  on  the  bed  and  bank  at  the  outer  bank  in  a 
meander.  But  in  practical  terms  the  boundary  shear  stress  is  a 
particularly  difficult  parameter  to  predict  accurately.  Indeed,  none 
specialists  even  find  it  difficult  to  visualize  boundary  shear  stress. 
Consequently,  it  is  desirable  to  relate  the  severity  of  bank  attack  and  toe 
scour  to  less  obscure  flow  descriptor,  such  as  near-bank  velocity.  Some 
modelers  even  prefer  to  relate  bank  attack  and  retreat  rates  to  near  bank 
velocities  instead  of  bank  shear  stress  (Odgaard,  1990).  Theory  shows 
that  near-bank  velocity  and  boundary  shear  stress  are  in  an,  case  closely 
related,  although  the  relation  between  them  is  neither  simple,  or  easily 
quantified  for  real  world  situations. 

The  preferred  treatment  to  stabilize  and  protect  the  outer  bank  in  a 
meander  bend  uses  a  blanket  of  loose  stone  called  riprap.  When  using 
riprap  it  is  necessary  to  select  the  appropriate  size  for  the  stone  on  the 
basis  of  the  intensity  of  flow  attack  as  represented  by  either  the  boundary 
shear  stress  on  the  outer  bank  or  the  flow  velocity  over  the  toe  of  the 
outer  bank.  Presently,  this  achieved  using  semi-empirical  diagrams  (Figs. 
1  and  2). 

The  first  (Fig.  1)  predicts  the  ratio  of  velocity  over  the  outer  bank 
toe  to  average  velocity  in  the  approach  channel  (Vtoe/Vavg)  as  a  function 
of  the  radius  of  curvature  to  width  ratio  for  the  bend  (Rc/w).  The  second 
(Fig.  2)  predicts  the  ratio  of  outer  bank  shear  stress  to  average  boundary 
shear  stress  in  the  approach  channel  (tb/to)  as  a  function  of  the  radius  cf 
curvature  to  width  ratio  for  the  bend. 
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Fig.  1  WES  design  diagram  for  prediction  of  outer  bank  velocity  at  a  bend 


The  velocity  diagram  uses  a  logarithmic  scale  for  the  independent 
variable  (Rc/w)  and  a  linear  scale  for  the  dependent  variable  (Vtoe/Vavg), 
Two  lines  are  plotted,  corresponding  to  natural  channels  (with 
asymmetrical  cross-sections)  and  trapezoidal  channels  (with  symmetrical 
cross-sections),  respectively.  The  ratio  of  outer  bank  to  mean  velocity  is 
markedly  higher  in  natural  than  trapezoidal  channels.  Plotted  as  straight 
lines  on  a  semi-log  graph,  these  lines  indicate  logarithmic  relations 
between  (Rc/w)  and  (Vtoe/Vavg)  for  the  two  types  of  channel.  The 
equations  of  the  lines  are  not  given,  but  analysis  of  the  graph  suggests 
that  they  approximate  to: 

Natural  Channels 

V55-  =  1.75  -0.5  log  (£)  (1) 


Trapezoidal  Channels 

-1.6  -0.71  log  (-) 
vavc  w 


(2) 


The  shear  stress  diagram  uses  logarithmic  axes  for  both  independent 
(Rc/w)  and  dependent  (tb/tQ)  variables.  Again,  two  lines  are  plotted,  this 

time  corresponding  to  smooth  and  rough  channels.  All  data  appear  to  come 
from  laboratory  flumes,  no  data  from  natural  rivers  are  included.  Rough 
channels  are  found  to  have  significantly  higher  stress  ratios  than  smooth 
channels,  for  the  same  value  of  (Rc/w),  although  the  line  for  rough 
channels  is  fitted  to  only  two  points  and  is  heavily  extrapolated.  Plotted 
as  straight  lines  on  log-log  graph,  these  lines  indicate  power  function 
relations  between  (Rc/w)  arvl  (tb/t0).  The  equation  for  the  smooth  channel 
line  is  given  on  the  diagram  as:- 


(3) 


No  equation  for  the  rough  channel  line  is  given,  but  examination  of  the 
graph  suggests  that  the  line  may  be  described  by:- 


t 

t 


b 

o 


(4) 


While  either  diagram  can  give  reasonable  results  when  used  with 
sound  engineering  judgement  and  with  careful  consideration  of  the  limits 


to  its  applicability,  it  is  nonetheless  desirable  to  develop  improved 
procedures  that  better  account  for  the  parameters  of  flow  hydraulics, 
boundary  roughness  and  channel  geometry  that  are  believed  to  influence 
flow  intensity  at  the  outer  bank  in  a  meander  bend.  Several  other  aspects 
of  bend  geometry,  channel  shape  and  boundary  roughness  have  been  shown 
to  influence  bend  flow  patterns  significantly  on  both  theoretical  and 
practical  grounds  (Thorne,  1978;  Hooke  and  Harvey,  1983;  Rais,  1984; 
Lapointe  and  Carson,  1986;  Pizzuto,  1987;  Thorne  and  Osman,  1988; 
Odgaard,  1989),  and  a  methood  which  uses  only  a  single  parameter  to 
characterize  the  bend,  ignoring  all  others,  cannot  account  for  ti  ase 
effects. 


Objectives 

The  objectives  of  this  study  are  to  develop  improved  analytical 
techniques  to  estimate  the  velocity  and  shear  stress  distributions  at  the 
outer  bank  in  a  river  bend.  The  approach  adopted  is  to  examine  these 
distributions  as  functions  of  the  planform  and  cross-sectional  geometry 
of  the  bend,  the  nature  of  the  bed  and  bank  materials,  and  the  planform  and 
average  flow  parameters  in  the  approach  channel. 

The  primary  objective  is  to  concentrate  on  defining  maximum  values 
of  depth  averaged  velocity  that  occur  in  the  bend  along  the  outer  bank  (that 
is  over  the  toe  of  revetted  banks).  The  second  aim  is  to  produce  the 
equivalent  relationships  for  boundary  shear  stress  at  the  outer  bank  in  a 
meander  bend. 

Emphasis  is  placed  on  basing  the  relationships  on  parameters  readily 
available  to  design  engineers,  rather  than  variables  such  as  "centerline 
mean  velocity’  which  although  theoretically  significant,  are  usually 
unkown  and  which  would  themselves  be  difficult  to  predict  or  estimate. 


Approaches  Adopted 

Broadly,  two  approaches  have  been  used.  The  first  is  based  op 
statistical  analysis  of  a  data  base  on  bend  flow  assembled  from  published 
and  unpublished  reports  of  studies  made  on  rivers  and  in  laboratory  flumes 
all  over  the  world.  The  second  attempts  a  more  theoretical  approach, 
being  based  on  application  of  three  recently  developed  mathematical 
models  of  bend  flow  hydraulics.  There  are  advantages  and  disadvantages 
to  both  approaches  and  these  are  discussed  in  the  sections  concerned  with 
the  Final  Discussion  and  Conclusions. 
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Data-Based  Approach 


Sources  of  Data 

Data  were  obtained  from  a  number  of  diverse  sources.  The  sources 
actually  used  are  listed  in  Appendix  A.  The  initial  data  came  from  studies 
undertaken  by  the  Principal  Investigators  and  their  colleagues  at  Colorado 
State  University,  London  University,  UK  and  the  University  of  East  Anglia, 
UK.  These  data  were  readily  to  hand  and  included  all  of  the  parameters 
necessary  for  this  analysis.  They  required  only  a  little  time  and  effort  to 
assemble. 

The  second  source  of  data  was  from  researchers  known  to  be  working 
on  bend  flow  problems  and  with  whom  the  Principal  Investigators  have 
good  working  relationships.  In  response  to  requests  from  the  Pi's  or  their 
research  associates,  copies  of  research  reports  and  published  articles 
containing  full  data  sets  were  supplied  by  these  individuals,  mostly  in  a 
timely  fashion.  This  allowed  easy  extraction  of  the  relevant  parameters. 
In  cases  where  a  particular  measurement  was  not  reported,  telephone 
calls  to  the  original  researchers  usually  elicited  the  missing  information. 

The  third  source  of  data  was  from  papers  published  in  professional 
and  learned  journals.  This  proved  to  be  the  least  satisfactory  source. 
Journal  papers  almost  never  contain  full  data  sets,  and  published  summary 
diagrams  of  the  distribution  of  parameters  such  as  depth-averaged 
velocity  are  too  small  to  be  used  for  data  extraction  with  any  degree  of 
accuracy  or  precision.  The  addresses  given  in  articles  are  often 
incomplete  or  out  of  date  and  telephone  and  FAX  numbers  are  omitted.  Most 
authors  were  extremely  slow  to  respond  to  written  enquiries  sent  by 
ordinary  mail  and  some  seemed  reluctant  to  part  with  data  at  all.  These 
problems  led  to  several  promising  leads  being  reluctantly  abandoned  and 
data  sets  excluded  from  the  analysis. 

The  data  set  which  has  resulted  is  then  not  universal  in  its  scope.  It 
does,  however,  contain  only  data  which  the  Principal  Investigators 
opinions  is  sound  and  complete.  The  range  of  sizes  and  types  of  channel 
encompassed  is  large  and  there  is  a  sufficient  number  of  entirely 
independent  data  sets  to  support  the  statistical  analysis.  Consequently,  it 
is  probable  that  the  addition  of  a  few  further  data  is  unlikely  to 
materially  alter  the  overall  distribution  of  data  or  the  outcome  of  the 
analyses. 

Data  base 

The  basic  data  assembled  in  this  study  are  listed  in  Tables  1,  2  and  3, 
for  Natural  Rivers,  Trapezoidal  Channels  and  Rectangular  Channels 
respectively.  The  published  and  unpublished  sources  of  data  are  listed 
separately  in  the  reference  section  of  this  report. 
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TABLt  I  •  BASIC  DATA  FOR  NATURAL  RIVERS 
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TABLE  J  -  BASIC  DATA  TOR  RECTANGULAR  CHANNELS 
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The  basic  data  were  used  to  derive  parameters  of  bend  geometry  and 
hydraulic  roughness  which  could  affect  the  pattern  of  flow  through  the 
bend.  The  derived  data  are  listed  in  Tables  4,  5  and  6  for  Natural  Rivers, 
Trapezoidal  Channels  and  Rectangular  Channels  respectively. 

An  important  aspect  of  any  experimentally  based  study  is  to  identify 
the  range  of  each  variable  observed.  When  applying  relationships  based  on 
the  experimental  results,  these  ranges  must  set  the  limits  to  the 
applicability  of  the  relations.  It  is  highly  speculative  and  very  risky  to 
apply  any  empirical  relationship  outside  the  range  of  data  from  which  it 
has  been  developed  and  tested.  The  range  of  each  of  the  variables  is  listed 
in  Data  Tables  7,  8  and  9  for  Natural,  Trapezoidal  and  Rectangular 
Channels,  respectively. 
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TABLE  4  -  DERIVED  DATA  FOR  NATURAL  RIVERS 
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TABLE  5  •  DERIVED  DATA  Kill  TRAPEZOIDAL  CHANNELS 
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Table  7  -  Range  of  Variables  for  Natural  Channels 


Measured  Variables 


Variable 

Units 

Range 

Radius  of  Curvature 

meters 

8  -  4,525 

Bend  Length 

meters 

16  -  5,633 

Width 

meters 

4  -  232 

Average  Depth 

meters 

0.4  -  5.65 

Outer  Bank  Angle 

degrees 

21  -  90 

Outer  Bank  Roughness 

Rough-Intermediate 

Median  Bed  Material  Size 

millimeters 

0.3  -  63 

Bedforms 

-  - 

Plane  -  Dunes 

Approach  Channel 

Straight-Meandering 

Average  Velocity 

meters/second 

0.42  -  1.47 

Depth-averaged  Toe  Velocity 

meters/second 

0.55  -  1.81 

Derived  Variables 


R/W 

0.75 

-  21.6 

L/w 

-  - 

1.45 

-  26.9 

w/d 

-  - 

9.05 

-  46.1 

d/D50 

-  - 

13.8 

-  18,833 

Vtoe/Vavg 

•  ” 

1.04 

-  1.57 

Table  8  -  Range  of  Variables  for  Trapezoidal  Channels 


Measured  Variables 

Variable 

Units 

Range 

Radius  of  Curvature 

meters 

1.5  -  15.24 

Bend  Length 

meters 

1.27  -  41.18 

Width 

meters 

0.61  -  6.76 

Average  Depth 

meters 

0.07  -  0.78 

Outer  Bank  Angle 

degrees 

27  -  56 

Outer  Bank  Roughness 

Smooth-Intermediate 

Median  Bed  Material  Size 

millimeters 

Smooth  -  38.1 

Bedforms 

Plane-Dunes 

Approach  Channel 

Straight-Meandering 

Average  Velocity 

meters/second 

0.34  -  1.07 

Depth-averaged  Toe  Velocity 

meters/second 

0.42  -  1.50 
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Derived  Variables 


R/w 

.  «. 

1.23  -  4.16 

L/w 

-  - 

1.31  -  16.88 

w/d 

-  - 

5.70  -  24.40 

d/D50 

11.0  -  2290 

Vtoe/Vavg 

-  - 

1.03  -  1.42 

Ttoe/Tavg 

-  “ 

1.00  -  2.80 

Table  9  -  Range  of  Variables  for  Rectangular  Channels 

Measured  Variables 


Variable 

Units 

Range 

Radius  of  Curvature 

meters 

0.8  -  50 

Bend  Length 

meters 

0.42  -  72.0 

Width 

meters 

0.30  -  6.00 

Average  Depth 

meters 

0.05  -  0.27 

Outer  Bank  Angle 

degrees 

18  -  90 

Outer  Bank  Roughness 

Rough-Smooth 

Median  Bed  Material  Size 

millimeters 

Smooth  -  40 

Bedforms 

Plane-Dunes 

Approach  Channel 

Straight 

Average  Velocity 

meters/second 

0.24  -  0.66 

Depth-averaged  Toe  Velocity 

meters/second 

0.35  -  0.81 

Derived  Variables 

R/w 

•  - 

0.83  -  9.01 

L/w 

— 

0.44  -  19.55 

w/d 

-  - 

2.22  -  24.0 

d/D50 

-  - 

4.3  -  2,500 

Vtoe/Vavg 

1.10  -  1.38 

Ttoe/Tavg 

1.20  -  2.46 

Examination  of  Data 

Before  undertaking  any  advanced  analysis  or  statistical  treatment  of 
data,  it  is  important  to  examine  the  data  carefully  in  the  light  of  existing 
knowledge  annd  theory.  This  allows  the  researcher  to  identify  expected 
and  unexpected  trends  and  relationships,  and  establishes  the  analytical 
framework  for  the  formal  treatment  of  the  data.  This,  fairly  lengthy, 
procedure  is  essential  if  the  resulting  relationships  are  to  have  physical 
as  well  as  statistical  significance. 
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The  first  step  was  to  establish  how  the  data  collected  in  this  study 
plotted  in  relation  to  the  design  curve  developed  by  the  US  Army  Engineer 
Waterways  Experiment  Station.  Hence,  a  semi-logarithmic  plot  of  (Rc/w) 
versus  (Vtoe/Vavg)  was  produced  for  the  Natural  River  data,  with  the  WES 
design  curve  marked  on  (Fig.  3a).  The  design  line  does  not  pass  trough  the 
points,  but  does  form  a  good  upper  bound  to  the  data  with  the  exception  of 
only  three  out  of  34  points.  Thus,  it  may  be  concluded  that  the  WES  design 
curve  represents  a  reasonable,  but  rather  conservative  approach  to  the 
estimation  of  (Vtoe/Vavg)  in  natural  channels.  This  is  essential  so  that  in 
the  final  design,  the  size  of  rip.ap  specified  is  always  on  the  safe  side.  A 
regression  line  through  the  scatter  of  the  points  for  Vtoe/Vavg  could  be 
used,  but  this  would  require  that  a  factor  of  safety  be  introduced  in  the 
relationship  between  the  critical  local  velocity  for  entrainment  and  the 
size  of  stone  used  in  a  revettment.  Present  WES  preference  is  to  position 
the  design  line  as  an  upper  bound  to  the  data,  so  that  all  of  the  zone  of 
uncertainty  is  on  one  side  of  the  line  (Oswald,  personal  communication, 
March  1990). 

However,  there  is  considerable  scatter  in  the  data,  and  this  deserves 
comment.  Partly,  it  is  a  result  of  the  methods  used  to  collect  the  data. 
Usually,  velocities  were  measured  at  a  finite  number  of  cross-sections 
around  each  bend.  In  some  studies  many  sections  were  used  (up  to  seven 
per  bend),  but  in  others  only  a  few  (less  than  three)  were  used.  Outer  bank 
velocities  at  intermediate  points  between  sections  were  not  measured. 
Consequently,  there  is  no  guarantee  that  the  actual  maximum  outer  bank  in 
a  bend  would  be  observed  in  any  study.  Indeed,  in  studies  with  only  a  few 
sections,  it  is  highly  probable  that  the  outer  bank  maximum  velocity  for  a 
bend  would  not  be  measured.  It  is  therefore  to  be  expected  that  field  data 
should  plot  either  close  to  or  below  a  line  defining  the  maximum  possible 
ratio  of  outer  bank  to  average  velocity.  However,  even  for  bends  with 
multiple  measured  sections,  the  data  often  plot  well  below  the  WES  line. 
This  suggests  that  there  may  be  further  variables  affecting  the  velocity 
ratio  which  are  unaccounted  for  in  the  WES  analysis. 

Points  for  bends  of  very  low  Rc/w  values  reveal  that  the  monolonic 
increase  in  Vtoe/  Vavg  observed  as  Rc/w  decreases  may  cease  at  an  Rc/w 
of  about  2.  For  Rc/w  values  less  than  2,  the  data  show  a  wide  range  of 
Vtoe/  Vavg  values,  but  the  velocity  ratio  never  exceeds  1.6.  This  accords 
with  other  recent  studies  of  bend  flow  in  very  tightly  curved  bends,  which 
has  shown  that  both  outer  bank  scour  pool  depth  and  outer  bank  retreat 
rate  may  actually  decrease  with  decreasing  Rc/w  for  bends  with  Rc/w 
less  than  2  (Biedenharn  et  al.,  1989;  Thorne,  1989).  This  is  not  unexpected 
theoretically,  as  there  is  a  major  discontinuity  in  the  way  the  pattern  of 
bend  flow  responds  to  increasing  bend  tightness  at  Rc/w  of  between  2  and 
3  (Bagnold,  1960).  Further  data  and  analyses  are  required  to  confirm  this 
tentative  finding. 
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It  Is  concluded  that  the  actual  ratio  of  outer  bank  toe  velocity  to 
average  velocity  at  a  bend  increases  as  the  ratio  of  radius  of  curvature  to 
width  decreases,  in  bends  with  Rc/w  greater  than  2.  In  a  natural  channel 
the  actual  velocity  ratio  observed  in  the  field  is  unlikely  to  exceed  the 
value  predicted  from  the  WES  design  curve,  but  it  is  likely  to  be 
considerably  lower  under  some  circumstances.  For  very  tight  bends  with 
Rc/w  less  than  2,  a  wide  range  of  Vtoe/Vavg  values  is  possible,  but 
maximum  values  never  exceed  1.6. 

Effect  of  Channel  Shape 

Figure  3b  shows  the  same  plot  for  trapezoidal  channels,  again  with 
the  relevant  WES  design  curve  superimposed.  The  trend  of  the  line  is 
clearly  correct,  but  the  data  tend  to  scatter  about  the  line  rather  than 
lying  near  or  below  it  as  in  the  case  of  natural  channels.  Three  out  of  ten 
points  lie  significantly  above  the  line,  suggesting  that  it  might  be  prone  to 
underestimating  the  actual  ratio  of  toe  to  average  velocity  under  some 
circumstances. 

Figure  3c  shows  the  same  plot  for  rectangular  channels.  Both  the 
lines  for  natural  and  trapezoidal  channels  are  superimposed.  The  data  tend 
to  plot  around  the  line  for  trapezoidal  channels,  eleven  points  lie  above 
and  six  below  the  line.  As  the  shape  of  a  rectangular  channel  is  something 
between  trapezoidal  and  natural,  the  plotting  position  of  the  points  is  as 
expected.  The  plot  suggests  that  Vtoe/Vavg  values  in  rectangular  channels 
are  lower  than  those  found  in  natural  channels,  but  may  be  somewhat 
higher  than  those  found  in  trapezoidal  channels. 

In  order  to  establish  which  other  variables  influence  the  velocity 
ratio  for  a  bend,  separate  semi-logarithmic  graphs  were  plotted  for 
further,  different  channel  characteristics. 

Effect  of  Bank  Roughness 

Figures  4a  and  4b  show  the  Rc/w  versus  Vtoe/Vavg  relations  for 
natural  bends  with  intermediate  roughness  outer  banks  and  rough  outer 
banks,  respectively.  Examination  of  the  plots  shows  complete  overlap 
between  the  data  clouds  for  the  two  bank  types.  This  suggests  that,  for  the 
range  of  bank  roughness  represented  in  the  bends  studied,  the  roughness  of 
the  outer  bank  did  not  materially  affect  the  velocity  ratio. 

The  banks  of  the  laboratory  flumes  used  to  generate  the  data  for 
trapezoidal  and  rectangular  channels  showed  an  insufficient  range  of 
roughness  to  allow  separation  of  the  data  in  this  way. 

Effect  of  Bedforms 

Figures  5a,  b  and  c  show  the  Rc/w  versus  Vtoe/Vavg  relations  for 
natural  bends  with  plane,  ripple  and  dune,  and  dune  bedforms,  respectively. 
Examination  of  the  plots  shows  complete  overlap  of  the  data  clouds  for  the 
three  bedforms,  suggesting  that  in  natural  Dends  the  bedform  did  not 
significantly  affect  the  velocity  ratio. 
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Fig.  5c  Natural  Rivers  -  Dune  Beds 
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Figures  6a  and  b  show  the  relations  for  rectangular  channels  with 
plane  and  dune  beds,  respectively.  The  plots  show  that  dune-bedded  bends 
seem  to  have  higher  velocity  ratios  than  plane-bedded  bends  for  the  same 
Rc/w  value.  This  appears  to  conflict  with  the  finding  for  natural  channels, 
that  bedform  did  not  affect  the  velocity  ratio. 

The  probable  explanation  centers  on  the  state  of  mobility  of  the  bed 
in  the  flume  studies.  The  presence  of  dunes  clearly  indicates  mobile-bed 
conditions  and,  therefore,  the  potential  for  scour  and  fill.  In  a  mobile-bed 
experiment  I  would  expect  a  scour  pool  to  develop  in  the  bed  adjacent  to 
the  outer  bank  and  a  point  bar  to  form  at  the  inner  bank,  so  that  in  time  the 
cross-section  would  come  to  some  extent  to  resemble  that  of  a  natural 
channel.  Data  collected  at  that  time  would  plot  between  the  lines  for 
fixed  and  natural  channels. 

Conversely,  a  plane  bed  indicates  immobile  conditions  and  a  fixed, 
rectangular  cross-section.  Data  from  such  a  channel  should  plot  close  to 
the  line  for  trapezoidal  channels.  Viewed  in  this  light,  the  data  for 
different  bedforms  simply  illustrate  the  effect  of  cross-sectional  shape, 
with  mobile-bed  rectangular  channels  approximating  to  natural  channels 
and  immobile-bed  rectangular  channels  being  similar  to  trapezoidal 
channels.  Re-examination  of  the  data  for  trapezoid?.1  channels  shows  that 
this  effect  is  also  evident  there.  The  two  points  which  plot  well  above  the 
WES  design  curve  both  come  from  a  channel  with  a  granular,  deformable 
bed. 

On  this  basis  it  seems  sensible  to  consider  laboratory  channels  as 
either  having  a  mobile  or  an  immobile  bed,  pooling  together  those  with 
initially  trapezoidal  and  rectangular  cross-sectons. 

Effect  of  Entrance  Conditions 

Figures  7a  and  b  show  the  Rc/w  versus  Vtoe/Vavg  relations  for 
natural  channel  bends  with  straight  and  meandering  entrance  conditions, 
respectively.  Examination  of  the  plots  suggests  that  for  the  same  Rc/w 
value,  bends  downstream  of  straight  reaches  have  higher  velocity  ratios 
than  those  downstream  of  meandering  reaches.  The  WES  design  curve 
forms  a  reasonable  upper  bound  to  the  data  for  bends  with  straight 
entrance  conditions,  but  significantly  over-estimates  the  increase  in  the 
velocity  ratio  that  accompanies  a  decrease  in  Rc/w  for  bends  in 
meandering  reaches.  The  discrepancy  increases  as  the  Rc/w  decreases. 

This  is,  potentially,  an  important  finding  because  it  suggests  that 
different  design  approaches  might  be  appropriate  for  single,  isolated 
bends  and  the  consectutive  bends  of  a  meandering  river. 

The  difference  may  arise  due  to  the  contrasting  transverse 
distributions  of  longstream  velocity  at  the  entrance  of  bends  with 
straight  and  meandering  reaches  upstream. 

At  the  end  of  a  long,  straight  approach  reach  the  maximum  velocity 
filament  is  close  to  the  channel  centerline.  It  must  cross  only  half  the 
channel  width  before  encountering  the  outer  bank  zone,  and  elevating  outer 
bank  toe  velocities  relative  to  the  average  velocity.  This  is  usually 
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Fig.  7a  Natural  Rivers  -  Straight  Entrance 
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achieved  just  after  the  bend  apex,  in  the  second  half  of  the  bend,  where 
skewing  of  the  flow  by  the  bend  is  at  its  strongest. 

Conversely,  at  the  exit  of  an  upstream  bend  the  maximum  velocity 
filament  is  located  adjacent  to  the  outer  bank  for  that  bend.  If  the  next 
bend  (which  is  of  opposite  curvature)  is  immediately  downstream,  then  at 
its  entrance  the  maximum  velocity  filament  is  located  near  the  inner  bank. 

It  must  cross  the  whole  channel  width  before  encountering  the  outer  bank 
zone.  This  requires  fully  developed  secondary  flow  and  a  long  bend,  and  is 
seldom  achieved  until  downstream  of  the  bend  exit,  where  the  strength  of 
skewed  flow  is  already  declining.  As  a  result,  the  ratio  of  outer  bank  to 
average  velocity  in  such  a  bend  may  be  lower  than  that  in  an  equivalent 
bend  at  the  end  of  a  straight  reach.  This  preliminary  finding  is  consistent 
with  long-held  ideas  on  the  effect  o *  entrance  conditions  on  bend  flow 
(see  for  example,  Chacinski  and  Francis,  1952)  and  it  merits  further 
research. 

The  only  break-down  of  the  data  to  show  a  systematic  and  significant 
difference  in  the  relationship  between  Rc/w  and  Vtoe/Vavg  was  that  by 
entrance  condition.  It  was,  therefore,  decided  to  continue  this  breakdown 
when  examining  the  impact  of  the  other  parameters  of  bend  geometry,  and 
bed  and  bank  characteristics. 

Effect  of  Bend  Length 

This  was  investigated  using  an  dimensionless  index  of  bend  length 
(L/w)  similar  in  concept  to  the  dimensionless  bend  curvature  index  Rc/w. 
Figures  8a  and  b  show  the  relations  for  natural  channel  bends  with 
straight  and  meandering  entrance  conditions,  respectively. 

The  trend  in  both  graphs  is  for  Vtoe/Vavg  to  increase  as  L/w,  the 
dimensionless  length  of  the  bend,  decreases.  This  seems  surprising  since, 
from  first  principles,  I  would  have  expected  that  the  asymmetry  of  the 
flow  would  be  greater  in  a  long  bend  than  a  short  one  due  to  the  greater 
length  ol  curved  channel  over  which  the  skewed  flow  may  develop. 
However,  it  must  be  noted  that  bend  length  is  intimately  related  to  the 
Rc/w  value  as  well  as  the  L/w  value. 

It  is  a  geometric  fact  that  tight  bends  must  be  short  in  length, 
otherwise  a  neck  cut-off  occurs.  Consequently,  there  is  a  high  positive 
correlation  (R  *  0.88)  between  Rc/w  and  L/w  (Table  10).  Hence,  short 
bends  have  high  Vtoe/Vavg  ratios  not  because  they  are  short,  but  because 
they  are  tight.  The  effects  of  bend  length  cannot  easily  be  separated  from 
those  of  radius  of  curvature  to  width  ratio.  Physically,  what  happens  is 
that  in  a  short,  tight  bend  the  maximum  velocity  filament  crosses  the 
channel  rather  abruptly  and  meets  the  outer  bank  at  an  acute  angle.  The 
result  is  that  acceleration  of  the  outer  bank  velocity  in  a  short,  tight  bend 
would  be  expected  to  be  greater  than  that  in  a  long/gentle  bend.  The  only 
way  to  separate  the  effect  of  L/w  acting  alone  is  through  multiple 
regression,  which  examines  the  variation  of  Vtoe/Vavg  caused  be  a  change 
in  L/w  after  the  effects  of  Rc/w  have  been  accounted  for.  This  rnalysis  is 
performed  later  in  the  report. 
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Table  10  Correlation  Matrix  Between  derived  Parameters 

Correlation  Matrix  for  Varlablea:  X-|  ...  X5 


Rc/w 
L/w 
w/d 
d/DSO 
Bank  Angle 


Rc/w  L/w _ w/d _ d/D50  BankAn... 


1 

.881 

1 

.758 

.688 

1 

.858 

.834 

.885 

1 

-.775 

-.777 

-.67 

-.86 

1 

The  results  for  trapezoidal  and  rectangular  channels  could  not  be 
broken  down  between  straight  and  meandering  entrance  conditions  due  to 
limitations  in  the  scope  of  the  data  set.  There  were  not  enough  cases  with 
meandering  entrance  conditions  to  allow  a  meaningful  analysis.  As  most 
of  the  data  come  from  straight  approach  conditions,  the  results  should  be 
on  the  safe  side  with  regard  to  meandering  approach  channels. 

The  results  are  plotted  in  Figs.  9  and  10  for  trapezoidal  and 
rectangular  channels,  respectively.  In  both  cases  the  trend  of  the  data  is 
apparently  similar  to  that  for  natural  channels  (Fig.  8a  and  b).  However, 
the  comparison  of  the  plots  with  those  for  natural  channels  emphasizes 
the  relatively  restricted  range  of  L/w  ratios  found  in  the  laboratory 
channels,  which  makes  it  less  easy  to  identify  any  trend  in  the  data. 

Effect  of  Aspect  Ratio 

Figures  11a  and  b  show  the  relations  for  aspect  ratio  (w/d)  versus 
Vtoe/Vavg  in  natural  channel  bends  with  straight  and  meandering  entrance 
conditions,  respectively. 

The  trend  in  both  graphs  is  for  Vtoe/Vavg  to  increase  as  width  to 
depth  ratio  decreases.  This  is  explained  by  the  different  flow  patterns  in 
narrow,  deep  channels  versus  shallow,  wide  channels.  In  a  shallow,  wide 
channel  bed  roughness  dominates  the  flow  pattern,  flow  is  pseudo 
two-dimensional,  and  both  curvature  and  outer  bank  effects  are  relatively 
less  important.  Conversely,  in  a  narrow,  deep  channel  the  flow  is  fully 
three-dimensional  and  strong  skew-induced  and  outer  bank  secondary 
currents  dominate  the  flow  pattern.  It  is  these  secondary  currents  which 
are  responsible  for  the  accelaration  of  near  bank  primary  velocities  in  the 
outer  bank  region.  Consequently,  the  increase  in  the  outer  bank  toe 
velocity  is  much  greater  in  a  narrow,  deep  channel  than  in  an  equivalent 
shallow,  wide  one.  Again,  Table  7  reveals  a  positive  correlation  between 
w/d  and  Rc/w  (R  *=  0.76),  however.  Therefore,  final  cosideration  of  the 
effect  of  w/d  on  Vtoe/Vavg  must  be  reserved  until  a  multiple  regression 
has  been  performed. 
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Fig.  11b  Natural  Rivers  -  Meandering  Approach  Channel 
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The  results  for  trapezoidal  and  rectangular  channels  are  shown 
in  Figs.  12  and  13,  respectively.  The  ranges  of  w/d  encountered  in 
laboratory  channels  are  severely  restricted  compared  to  natural  channels, 
but  the  trend  in  the  data  does  not  contradict  the  results  for  natural 
channels. 

Effect  of  Relative  Depth 

Figures  14a  and  b  show  the  relations  for  relative  depth  (d/D50) 
versus  Vtoe/Vavg  in  natural  channel  bends  with  straight  and  meandering 
entrance  conditions,  respectively. 

The  trend  in  both  graphs  is  for  Vtoe/Vavg  to  increase  as  relative 
depth  decreases.  This  is  explained  by  the  effect  of  flow  resistance  on 
average  velocity  in  a  channel.  Relative  depth  (the  inverse  of  relative 
roughness  (D50/d))  characterizes  how  deep  the  flow  is  in  terms  of  the 
median  size  of  the  bed  material.  It  represents  the  hydraulic  smoothness 
of  the  channel,  which  is  primarily  responsible  for  determining  the  flow 
resistance  coefficient,  and  hence  the  average  velocity.  In  a  relatively 
deep,  smooth  channel  (high  d/D50)  resistance  is  low  and  average  velocity 
high  when  compared  to  an  otherwise  equivalent  shallow,  rough  channel 
(low  d/D50).  Consequently,  when  the  outer  bank  toe  velocity  is  divided  by 
that  average  velocity  in  a  deep,  smooth  channel,  the  quotient  is  lower  than 
in  an  otherwise  equivalent  shallow,  rough  channel.  The  difference 
becomes  more  marked  as  d/D50  decreases.  It  should  be  noted  though  that 
this  effect  operates  on  the  quotient  Vtoe/Vavg.  It  does  not  mean  that 
absolute  values  of  Vtoe  will  be  lower  in  a  deep,  smooth  channel  than  in  a 
shallow,  rough  one,  because  the  actual  value  of  Vavg  may  also  be  much 
higher  in  the  smoother  channel.  The  correlation  matrix  (Table  7)  shows  a 
strong,  positive  correlation  between  d/D50  and  Rc/w  (R  •=  0.86).  Hence, 
final  judgement  on  the  relationship  bewteen  d/D50  and  Vtoe/Vavg  must 
await  the  multiple  regression  analysis. 

There  is  considerable  scatter  in  both  the  graphs,  especially  that  for 
straight  approach  conditions  (Fig.  14a),  but  the  points  for  bends  with 
straight  entrance  approach  channels  consistently  plot  higher  than  those 
for  meandering  approaches.  This  suggests  that  compared  to  meandering 
rivers,  maximum  values  of  Vtoe/Vavg  are  higher  for  bends  with  straight 
approach  channels,  and  increase  more  markedly  as  d/D50  decreases. 

The  results  for  trapezoidal  and  rectangular  channels  are  shown  in 
Figs.  15  and  16,  respectively.  The  trends  in  the  data  support  the  results 
for  natural  channels  over  a  wide  range  of  relative  depths. 

Effect  of  Outer  Bank  Angle 

Figures  17a  and  b  show  the  relations  for  outer  bank  angle  versus 
Vtoe/Vavg  in  natural  channel  bends  with  straight  and  meandering  entrance 
conditions,  respectively. 

The  trend  in  both  graphs  is  for  Vtoe/Vavg  to  increase  with  increasing 
bank  angle. 
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Fig.  14b  Natural  Rivera  -  Meandering  Approach  Channel 
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Fig.  1 7a  Natural  Rivers  -  Straight  Approach  Channel 
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Fig.  1 7b  Natural  Rivers  *  Meandering  Approach  Channel 
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There  is  weak  evidence  that  the  Vtoe/Vavg  ratio  might  level  off,  or 
even  decrease  again  as  bank  angle  increases  from  about  75  to  90  degrees. 
The  increase  in  Vtoe/Vavg  with  increasing  bank  angle  is  explained  by  the 
effect  of  outer  bank  angle  on  flow  patterns  adjacent  to  the  bank.  In  a 
channel  with  a  shelving  outer  bank,  depth  decreases  gradually  as  the 
water's  edge  is  approached.  Relative  roughness  (D50/d)  increases  to  high 
values,  and  momentum  (both  longstream  and  cross-stream)  is  lost  to 
friction  at  the  boundary.  The  high  roughness,  low  depth  and  low  velocities 
combine  to  suppress  all  but  the  stongest  secondary  currents,  and  no  outer 
bank  cell  of  secondary  circulation  develops.  Hence,  the  maximum  primary 
velocity  in  the  bank  zone  is  relatively  low,  and  it  stays  at  the  free 
surface,  well  alway  from  the  toe. 

But  in  a  bank  with  a  steep  outer  bank,  the  depth  increases  quickly  at 
the  water's  edge.  High  velocity  flow  is  able  to  move  in  close  to  the  bank, 
and  secondary  currents  are  pronounced.  The  junction  of  the  water  surface 
and  the  bank  acts  as  a  stagnation  point  where  the  velocity  of  outward  flow 
near  the  surface  must  go  to  zero  very  abruptly  as  the  skew-induced  cell 
meets  the  outer  bank.  This  results  in  a  small  cell  of  reverse  rotation. 
Combination  of  the  outer  bank  and  skew-incduced  cells  with  the  high 
velocity  primary  flow,  produces  a  depressed  maximum  in  the  region  above 
the  bank  toe,  which  can  produce  much  higher  Vtoe/Vavg  values  than  those 
found  adjacent  to  a  shelving,  low-angle  bank. 

The  apparent  reduction  in  Vtoe/Vavg  observed  in  bends  with  near 
vertical  outer  banks  probably  results  from  the  definition  of  Vtoe  used 
here.  This  is  the  depth-averaged  longstream  velocity  over  the  toe  of  the 
bank.  But,  when  the  bank  angle  is  close  to  90°,  this  definition  breaks  down 
because  the  bank  toe  and  bank  top  practically  coincide.  That  is,  the  Vtoe 
corresponds  to  the  velocity  at  the  water's  edge,  which  is  relatively 
reduced  by  simple  boundary  friction  compared  to  that  a  finite  distance  out 
into  the  channel.  Also,  the  maximum  near  bank  velocity  is  located  about 
one  third  of  the  way  up  the  bank,  rather  than  at  the  toe.  In  such  cases,  a 
modified  parameter  is  needed  to  represent  the  characteristic  velocity 
adjacent  to  the  outer  bank. 

The  correlation  between  bank  angle  and  Rc/w  (R  «  -0.78)  is  a  strong, 
negative  one  (Table  10).  This  indicates  that  steeper  banks  are  encountered 
in  tighter  bends.  Therefore  when  considering  the  effect  of  bank  angle  on 
Vtoe/Vavg,  final  conclusions  cannot  be  drawn  until  after  completion  of  the 
multiple  regression  analysis. 

The  scatter  of  points  for  straight  versus  meandering  entrance 
conditions  shows  that  the  distributions  overlap  for  low  bank  angles,  but 
that  the  increase  in  Vtoe/Vavg  with  increased  bank  angle  appears  more 
marked  in  channels  with  a  straight  approach. 

The  data  trapezoidal  channel  are  shown  in  Fig.  18.  Although  sparce, 
they  do  not  contradict  the  results  for  natural  channels.  Rectangular 
channels  have  vertical  outer  banks  by  definition  and  so  no  plot  could  be 
generated  for  this  variable. 
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Main  Points 


After  this  preliminary  examination  of  the  data  the  following  points 
were  noted: 

1.  The  WES  design  curve  represents  a  conservative  approach  to  the 
estimation  of  (Vtoe/Vavg)  in  natural  channels. 

2.  Points  for  bends  of  very  low  Rc/w  values  reveal  that  the  monotonic 
increase  in  Vtoe/  Vavg  observed  as  Rc/w  decreases  may  cease  at  an 
Rc/w  of  about  2. 

3.  For  trapezoidal  channels,  the  WES  design  curve  might  be  prone  to 
underestimating  the  actual  ratio  of  toe  to  average  velocity  under  some 
circumstances. 

4.  Vtoe/  Vavg  values  in  immobile  bed  rectangular  channels  are  similar  to 
those  found  in  immobile  bed  trapezoidal  channels. 

5.  The  roughness  of  the  outer  bank  does  not  significantly  affect  the  ratio 
of  outer  bank  toe  velocity  to  average  velocity  in  the  bends  studied. 

6.  The  presence  of  bedforms  does  not  significantly  affect  the  ratio  of 
outer  bank  toe  velocity  to  average  velocity  in  the  bends  studied.  It  is 
noted,  though,  that  the  velocity  ratio  for  immobile,  trapezoidal  or 
rectangular  channels  is  significantly  lower  than  that  for  mobile-bed, 
trapezoidal  or  rectangular  channels  with  the  same  Rc/w  ratio. 

7.  The  ratio  of  outer  bank  toe  velocity  to  average  velocity  in  a  bend 
immediately  downstream  of  a  straight  reach  may  be  significantly 
higher  than  that  in  a  bend  downstream  of  a  bend  of  opposite  curvature. 
The  WES  design  line  appears  to  be  a  good  upper  bound  for  bends 
downstream  of  straight  reaches  but  may  over-estimate  the  velocity 
ratio  in  bends  in  meandering  reaches. 

8.  Bend  geometry  parameters  L/w,  w/d  and  d/D50  are  all  show  strong, 
positive  correlations  with  Rc/w.  This  makes  it  difficult  to  identify 
their  individual  effects  on  Vtoe/Vavg.  From  the  individual  plots  it 
appears  that  Vtoe/Vavg  increases  as  L/w,  w/d  and  d/D50  decrease. 

9.  Bank  angle,  a,  shows  a  strong,  negative  correlation  with  Rc/w. 
Vtoe/Vavg  appears  to  increase  as  the  bank  angle  increases.  There  is 
some  evidence  that  Vtoe/Vavg  might  decrease  as  bank  angle  increases 
from  about  75  to  90  degrees. 

It  should  be  noted  that  there  is  considerable  scatter  in  the  graphs 
produced  here.  It  would  be  a  simple  matter  to  ’massage"  the  data  in  order 
to  obtain  tighter  distributions.  This  temptation  has  been  resisted, 
however,  and  the  data  are  exactly  as  they  were  extracted  from  the  various 
reports  and  papers.  No  doubt,  some  errors  of  judgement  have  been  made 
and  some  bends  have  been  mis-classified.  But  re-assigning  bends  to 
different  "straight"  and  "meandering"  categories  now  would  undermine  the 
objectivity  of  the  study.  The  analysis  is  based  on  the  our  best  judgement 
without  the  benefit  of  hind-sight.  If  it  cannot  produce  good  results  under 
these  circumstances,  then  it  would  be  unlikely  to  be  of  practical  use. 
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Analysis  of  Data 

The  examination  of  the  data  has  established  that  variables  other  than 
Rc/w  appear  to  influence  the  ratio  of  outer  bank  toe  velocity  to  average 
velocity  in  a  channel  bend.  Up  to  now,  each  of  the  potentially  significant 
variables  has  been  considered  separately,  as  if  it  acted  alone  to  influence 
Vtoe/Vavg.  In  nature  this  is  not  the  case,  as  the  actual  Vtoe/Vavg  value 
for  a  bend  is  the  product  of  the  mutual  interaction  of  all  the  significant 
controlling  variables.  Significant  correlations  exist  between  the 

variables,  and  these  may  obscure  the  true  effect  of  each  on  the  dependent 
variable  (Vtoe/Vavg). 

To  develop  a  predictive  approach  which  reflects  this  situation, 
multiple  regression  was  used  to  produce  equations  expressing  the 
dependent  variable  Vtoe/Vavg  as  a  function  of  all  "independent"  variables 
identified  as  being  potentially  significant  in  the  previous  section.  A 
stepwise  regression  was  used.  This  adds  each  independent  variable  in  turn 
to  the  regression  analysis,  rejecting  any  variable  which  does  not  add 
significantly  to  the  strength  of  the  regression  equation. 

From  the  examination  of  the  data  it  was  concluded  that  the  linear 
relation  between  Rc/w  and  Vtoe/Vavg  may  break  down  for  Rc/w  values 
less  than  2.  There  are  sound  theoretical  reasons  to  expect  this,  and  the 
data  are  consistent  with  the  idea  that  as  bends  tighten  to  Rc/w  values 
less  than  2  major  changes  in  flow  pattern  occur,  often  leading  to 

impinging  flow  and  areas  of  separation  at  the  outer  bank.  In  view  of  this, 
it  was  decided  to  limit  the  curve  fitting  to  bends  with  Rc/w  values  equal 
to  or  greater  than  2. 

It  was  noted  earlier  that  in  some  studies  only  a  few  sections  were 

monitored  in  each  bend,  and  that  consequently  the  data  collected  do  not 

represent  the  absolute  maximum  ratio  of  Vtoe/Vavg  for  that  bend.  This 
explains  why  some  points  plot  low  in  the  various  distributions  shown  in 
Fig.  3.  Taking  this  fact  together  with  the  recommendation  from  WES  that 
an  upper  bound  line  is  preferable  to  a  best-fit  line  when  predicting 
Vtoe/Vavg  for  riprap  design,  it  was  decided  to  perform  the  multiple 

regression  using  only  points  from  the  top  edge  of  the  scatter.  These 
turned  out  to  come  mostly  from  bends  which  were  intensively  studied  (for 
example,  studies  by  Thorne  et  al.,  Dietrich,  Bridge).  Some  data  came  from 
single  sections  too,  apparently  where  those  sections  happened  to  have 
coincided  with  the  highest  values  of  Vtoe/Vavg  for  the  bend.  Hence,  the 
correlation  coeffiecients  indicate  the  linearity  of  the  upper  surface  of  the 
data  cloud,  rather  than  the  strength  of  the  regression  as  such. 

Prediction  of  Outer  Bank  Velocity  in  Open  Channel  Bendways 

Natural  Rivers 

In  line  with  the  finding  that  the  approach  condition  did  appear  to 
affect  the  velocity  ratio  for  natural  rivers,  separate  analyses  were 
performed  for  straight  and  meandering  approach  conditions. 
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Natural  Rivers:  Straight  Approach  Conditions 

The  points  used  to  define  the  upper  boundary  are  indicated  in  Fig.  19. 
The  equation  for  natural  rivers  with  straight  entrance  conditions  is: 


TT5-  *  129  ■ 0  5  >°g  ©  *  0  4  to*  ©  -  0.06  log  (A)  +  0.16  log  (a)  (5) 

AVG  w  w  ux, 

The  adjusted  coefficient  of  determination  for  this  equation  (r2)  is 
0.905,  indicating  that  the  equation  is  well  fitted  to  the  upper  surface  of 
the  data  cloud.  The  width-depth  ratio  did  not  contribute  sufficiently  to 
appear  in  the  equation. 

For  comparison  with  the  WES  design  curve,  and  to  evaluate  whether 
the  improvement  in  accuracy  merits  the  increased  complexity  of  eqn.  5,  a 
line  was  also  fitted  to  the  data  using  simple  regression  for  Vtoe/Vavg  as 
a  function  of  Rc/w.  The  resulting  equation  is: 

y"~  =  166  * 042  J°g  (~)  (6) 

VAVG  w 

The  adjusted  coefficient  of  determination  for  this  equation  (r2)  is 
0.90,  indicating  good  linearity  at  the  upper  edge  of  the  data  cloud.  The 
resulting  line  is  shown  in  Fig.  19. 

Natural  Channels:  Meandering  Approach  Conditions 

The  points  used  to  define  the  upper  boundary  are  shown  in  Fig.  20. 
The  resulting  equation  for  natural  rivers  with  meandering  entrance 
conditions  is: 

V“^  *  095 '  °- 16  loS  (“)  -  0.06  log  (^)  +  0.06  log  (~)  +  0.25  log  (a )  (7) 

The  adjusted  coefficient  of  determination  for  this  equation  (r2)  is 
0.985,  indicating  that  the  equation  is  well  fitted  to  the  upper  surface  of 
the  data  cloud.  The  relative  depth  did  not  contribute  sufficiently  to  appear 
in  the  equation. 

For  comparison  with  the  WES  design  curve,  and  to  evaluate  whether 
the  improvement  in  accuracy  merits  the  increased  complexity  of  eqn.  7,  a 
line  was  also  fitted  to  the  data  using  simple  regression  for  Vtoe/Vavg  as 
a  function  of  Rc/w.  The  resulting  equation  is: 

^  -  1-4 -0.24 log©  (8) 
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The  adjusted  coefficient  of  determination  for  this  equation  (r2)  is 
0.87,  again  indicating  good  linearity  at  the  upper  edge  of  the  data  cloud. 
The  line  produced  by  this  equation  is  shown  in  Fig.  20. 

Trapezoidal  and  Rectangular  Channels 

There  was  no  obvious  difference  between  the  data  clouds  for 
trapezoidal  and  rectangular  channels,  but  the  state  of  bed  mobility  did 
appear  to  affect  the  outer  bank  to  average  velocity  ratio  for  laboratory 
channels.  The  plots  for  mobile  and  fixed  bed  channels  are  shown  in  Figs. 
21  and  22  respectively.  Hence,  separate  analyses  were  carried  out  for 
mobile  and  fixed-bed  channels. 

Laboratory  Channels  with  Mobile  Beds 

The  points  used  to  define  the  upper  boundary  are  indicated  in  Fig.  21. 
There  was  insufficient  range  in  the  other  variables  (L/w,  w/d,  d/D50,  a)  to 
allow  their  meaningful  inclusion  in  the  analysis.  Hence  a  predictive 
equation  based  on  Rc/w  was  produced.  This  is: 

=  1.55  *  0.41  log  <*> 

AVG 

The  adjusted  coefficient  of  determination  for  this  equation  (r2)  is 
0.98,  indicating  good  linearity  at  the  upper  edge  of  the  data  cloud.  The  line 
produced  by  this  equation  is  shown  in  Fig.  21. 

Laboratory  Channels  with  Immobile  Beds 

The  points  used  to  define  the  upper  boundary  are  indicated  in  Fig.  22. 
The  WES  diagram  (Fig.  1),  indicates  that  Vtoe/Vavg  is  independent  of  Rc/w 
for  Rc/w  values  greater  than  6.  Accordingly,  only  points  with  Rc/w  values 
less  than  6  were  used  to  define  the  upper  boundary.  There  was 
insufficient  range  in  the  other  variables  (L/w,  w/d,  d/D50,  a)  to  allow 
their  meaningful  inclusion  in  the  analysis.  Hence  a  predictive  equation 
based  solely  on  Rc/w  was  produced.  This  is: 

1.73  -  0.89  log®  (10) 

VAVG  W 

The  line  produced  by  this  equation  is  shown  in  Fig.  22. 
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Prediction  of  Outer  Bank  Shear  Stress  in  Open  Channel  Bendways 

Very  little  reliable  data  on  shear  stress  distributions  in  natural 
channels  could  be  obtained.  Consequently,  the  analysis  was  limited  to  data 
from  laboratory  flumes.  It  should  be  noted  that  in  very  few  cases  were 
both  velocity  and  shear  stress  data  available  from  the  same  series  of 
experiments.  Therefore,  the  data  points  used  in  this  part  of  the  study  do 
not  correspond  to  those  used  to  investigate  Vtoe/Vavg. 

Nevertheless,  experience  gained  in  the  analysis  of  Vtoe/Vavg  did 
suggest  that  the  state  of  bed  mobility  might  affect  the  shear  stress  ratio 
for  a  bend.  The  plots  for  mobile  and  fixed  bed  channels  are  shown  in  Figs. 
23  and  24  respectively.  Separate  analyses  were  performed  for  channels 
with  mobile  and  immobile  beds. 

There  was  insufficient  data  to  support  multiple  regression,  and  so 
analysis  was  limited  to  simple  regression  on  upper  boundary  points  to 
predict  tb/to  as  afunction  of  Rc/w. 

Laboratory  Channels  with  Mobile  Beds 

The  points  used  to  fit  a  line  to  the  data  are  indicated  in  Fig.  23.  The 
equation  of  the  resulting  line  is: 

t.  .n  ^ 

-*-3.3i£)  an 

t  'wy 

O 

which  is  almost  identical  to  the  equation  of  the  WES  design  curve  for 
rough  channels.  However,  this  line  is  based  on  very  limited  data.  As  a 
result  the  top  of  the  data  cloud  is  poorly  defined,  producing  an 
unacceptably  low  regression  coefficient  of  only  0.27  for  the  upper 
boundary  line.  The  line  is  shown  in  Fig.  23. 

Further  data  from  natural  and  mobile-bed  laboratory  channels  are 
urgently  needed  to  confirm  the  form,  exponent  and  constant  in  this 
equation. 

Laboratory  Channels  with  Immobile  Beds 

The  points  used  to  fit  a  line  to  the  data  are  indicated  in  Fig.  23.  The 
WES  diagram  limits  analysis  to  bends  with  Rc/w  values  between  1  and  5. 
For  consistency,  these  limits  were  applied  to  the  data  used  here  also.  The 
equation  of  the  resulting  line  is: 

t  b  -0.5 

— *  «  2.92  (— )  (12) 

t  vw 

o 

which  is  very  similar  to  the  WES  design  equation  for  smooth  channels.  The 
high  correlation  coefficient  r  «=  0.96  indicates  good  linearity  in  the  upper 
boundary  of  the  data. 
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Fig.  23  Rectangular  Channels  with  Mobile  Beds 


Modeling  Approach 


Background 

During  the  last  decade,  a  large  number  of  papers  presenting 
mathematical  models  for  predicting  the  flow  and  bed  topography  in  river 
bends  has  appeared  in  the  literature.  The  authors  have  expressed  a  variety 
of  opinions  regarding  the  relative  importance  of  the  factors  which  control 
meander  morphology.  Most  models  are  intended  for  use  by 
sedimentologists  interested  in  reconstructing  past  environments  (eg. 
Allen,  1970,  Bridge,  1976),  or  by  river  engineers  attempting  to  predict  the 
distribution  of  scour  as  an  aid  in  the  design  of  successful  channel 
stabilization  schemes  (eg  Odgaard,  1987).  However,  many  models  may  be 
criticized  because  they  are  unnecessarily  complicated  and  esoteric  for 
these  tasks,  and  because  no  attempt  is  made  by  the  author(s)  to 
recommend  where,  when,  and  now  they  should  be  applied  to  natural 
waterways.  Some  very  pertinent  remarks  were  made  concerning  the  role 
of  mathematical  models  in  fluvial  geomorphology  and  river  engineering 
during  a  discussion  by  conference  participants  at  the  concluding  session  of 
the  ASCE  Rivers  ’83  conference,  New  Orleans,  1983  (Elliott,  1983).  These 
are  encapsulated  by  a  comment  from  Charles  Neill,  a  practising  river 
engineer  from  Northwest  Hydraulic  Consultants,  Canada.  He  said, 

"It  is  important  that  mathematical  models  should  have  a  good 
familiarity  with  the  range  of  features  encountered.  ...It  would 
be  a  service  to  the  profession  if  these  could  be  used  to 
produce  generalized  tables,  graphs  or  programs  that 
would  enable  reasonable  estimates  of  velocity  and  shear 
distributions  to  be  made  by  practising  engineers,  without 
the  necessity  of  access  to  an  elaborate  modelling  facility"  (Neill, 
1983) 

Difficulties  are  often  encountered  when  attempting  to  apply  models 
to  natural  systems.  If  the  model  requires  certain  input  parameters  that 
need  to  be  known  or  to  be  measured  in  the  field  (the  centreline  mean 
velocity  or  the  mean  Darcy-Weisbach  friction  factor  are  examples),  it  is 
often  difficult  to  assign  a  value  with  confidence.  Often,  field 
measurements  unavailable  and  estimates  are  unreliable,  and  it  may  be  the 
case  that  the  model  output  may  highly  sensitive  to  incorrect  values  having 
been  assigned  to  the  input  parameters. 

Also,  many  models  are  written  by,  and  apparently  for,  researchers. 
Unless  you  are  a  specialist,  expert  in  programming,  three  dimensional 
fluid  mechanics,  and  mathematics,  it  is  virtually  impossible  to  use  them 
for  a  real  world  application  without  detailed  assisstance  from  the  author 
of  the  computer  coding. 

In  this  part  of  the  project,  we  examined  a  number  of  mathematical 
models,  hoping  to  select  several  to  try  as  predictors  of  outerbank  velocity 
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in  bendways.  Before  going  on  to  report  the  results  from  the  models 
selected,  it  is  relevant  to  present  a  short  outline  of  the  physical  basis  for 
the  models. 

Basic  Principles  of  Numerical  Modeling  of  Bend  Flow 

The  numerical  modelling  of  flow  and  sediment  processes  in  river 
bends  is  a  subject  that  is  receiving  ever-increasing  attention,  and  a  large 
number  of  models  are  available.  There  are  two  components  to  most 
models.  The  first  involves  a  solution  to  the  equations  of  motion  for  fluid 
flow.  The  second  is  the  interaction  between  the  flow  and  the  bed 
topography.  This  requires  balancing  the  different  forces  acting  on 
bed-material  particles  to  produce  an  equilibrium  bed  topography. 

All  flow  models  start  with  the  equations  of  motion  for  fluid  flow. 
For  application  to  bend  flow,  the  equations  are  usually  written  in 
cylindrical  coordinates.  Given  below  are  the  equations  of  motion  for  the 
steady  flow  of  an  incompressible  fluid  in  an  orthogonal  cartesian 
coordinate  system  (Rozovskii,  1957).  The  velocity  components  in  the  s 
(streamwise),  n  (perpendicular  to  s-axis),  and  z  (vertical  upwards  from 
the  stream  bed)  directions  are  denoted  u,  v,  and  w  respectively,  r  -  local 
radius  of  curvature;  p  «  pressure;  and  F  ■  friction  term  in  the  s,  n  and  z 
directions  respectively. 

du  du  3u  uv  1  dp  „ 

u —  +  v —  +  w- — b  —  - +  F. 

ds  dn  dz  r  r  ds 

dv  dv  dv  u2  1  dp  _ 

u —  +  v —  +  w - -  - - +  F 

ds  dn  dz  r  r  dn  n 

dw  dw  dz  1  dp  _ 

u —  +  v —  +  w —  +  g  - -  +  Fi 

ds  dn  dw  r  dz 

The  left  hand  sides  of  these  equations  are  the  convective  acceleration 
terms.  There  are  no  local  acceleration  terms  (d/dt)  and  so  strictly 
speaking  the  models  are  not  prognostic  but  diagnostic. 

The  continuity  equation  for  3-dimensional,  incompressible  flow  is 
also  specified: 

du  1  d(vr)  dw 

—  + - + =  0 

ds  r  dn  dz 

The  principal  cross-stream  and  downstream  force  balances  are  between: 

1)  the  centrifugal  and  pressure  gradient  forces  in  the  cross-stream 
direction:  and 

2)  the  downstream  balance  between  gravitational  and  frictional  forces. 
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Secondary  circulation  is  usually  considered  to  be  the  most  important 
effect  of  curvature  on  flow,  but  the  tilting  of  the  water  surface  is  also 
very  important,  because  it  alters  the  downstream  slope  of  the  water 
surface,  generating  large  cross-stream  variation  in  the  downstream 
boundary  shear  stress  and  velocity  fields. 

Flow  and  bed  topography  models  attempt  to  simulate  the  bed 
morphology  of  a  channel  bend  by  assuming  that,  at  equilibrium,  the  forces 
directed  inwards  and  outwards  on  each  bed  particle  are  balanced  (Allen, 
1970;  Bridge, 1977).  This  means  that  particles  of  different  sizes  travel 
along  paths  of  equal  depth  along  the  channel,  under  the  influence  of 
longstream  drag.  Models  that  use  this  scheme  differ  in  the  way  in  which 
the  forces  of  lift  and  longstream  drag  are  determined.  The  balance  of 
forces  acting  on  a  particle  in  the  mean  flow  direction  is: 

Fd  cos  d  =  (W-Fl)  cosa  tanj 

where  FD  ■  drag  force,  d  *=  deviation  of  the  bed  shear  stress  vector  from 
the  longstream  direction,  W  -  submerged  weight,  FL  «=  lift  force,  tanj  is 
the  dynamic  friction  coefficient  due  to  collisions  with  the  bed  and  other 
grains,  and  a  «  transverse  slope  of  the  point  bar  surface  (Bridge,  1977) 
(see  figure  25). 


DOWNSTREAM  FUJW 
w*ler  soface  DIRECTION 


Figure  25  Definition  diagram  for  flow  in  an  open  channel  bend  (adapted 
from  Bridge,  1977) 

The  transverse  force  balance  is  therefore: 

Fq  sin  d  =  (W-Fl)  sin  a 

For  any  given  point  on  the  transverse  bed  profile,  the  balance  of  drag  and 
immersed  weight  components  acting  on  a  bed  particle  is: 
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2  4  3 

p  (D/2)  ^  tan  d  =  —  p  (d/2)  '  (s  -  r)  g  sin  a 

Where  D  -  particle  diameter,  tx  ■=  longstream  bed  shear  stress,  and  s  and  r 
«  the  sediment  and  fluid  densities,  respectively. 

The  theory  assumes  that  the  particles  are  moving  as  contact  load.  It 
is  also  important  to  note  that  for  suspended  particles,  forces  due  to  lift 
and  the  cross-stream  component  of  particle  weight  are  insignificant. 
Another  assumption  common  to  many  models  is  that  the  angular  deviation 
(d)  of  the  shear  stress  vector  from  the  downstream  direction  (and 
therefore  the  local  transverse  bed  slope)  is  proportional  to  the  ratio 
between  the  depth  and  radius  of  curvature: 

.  h 
land- C  — 


where  C  is  an  empirical  coefficient.  This  relation,  developed  by  Rozovskii 
(1957)  from  the  equations  of  motion,  is  actually  applicable  only  to 
fully-developed  secondary  flow.  Fully-developed  flow  occurs  in  the 
downstream  part  of  long,  constant-radius  reaches  where  flow  and  bed 
topography  remain  constant  with  distance  downstream  and  are  independent 
of  upstream  conditions.  For  developing  flow  where  flow  and  bed 
topography  do  not  remain  constant  with  distance,  the  governing  equations 
are  more  complex  and  difficult  to  solve. 

As  the  outward  component  of  gravity  is  proportional  to  the  cube  of 
the  diameter  of  the  grain  whilst  the  inward-acting  drag  on  the  particle  is 
proportional  to  the  square  of  the  diameter.  This  leads  to  a  sorting 
mechanism,  recognized  as  an  important  process  in  meander  bends  by  Allen 
(1970),  Bridge  (1977),  Dietrich  and  Smith  (1984)  and  Parker  and  Andrews 
(1985),  whereby  for  the  same  velocity,  larger  particles  will  tend  to  roll, 
due  to  gravity,  out  towards  the  pool,  while  smaller  ones  will  tend  to  be 
swept  inwards  by  fluid  drag  (Fig.  26).  Wilson  (1973)  used  the  same 
principle  to  propose  a-*  explanation  for  sorting  in  straight  channels. 

Review  of  Important  Bendflow  Models 

Most  researchers  begin  their  analyses  with  the  equations  of  motion 
and  continuity,  then  simplify  them  until  a  solution  is  possible  (eg.  Dietrich 
and  Smith,  1983;  Odgaard,  1987).  Others  use  simpler  concepts  of 
momentum  and  force,  and  treat  the  flow  as  a  single  unit,  using 
depth-averaged  equations  (eg.  Dietrich,  1988).  The  problem  with  the 
latter  approach  is  that  the  effect  of  the  secondary  currents  cannot  be 
accounted  for  by  the  governing  equations. 
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Figure  26  Forces  acting  on  a  particle  at  or  near  the  bed.  Opposing 
frictional  resistance  of  the  bed  not  shown.  Forces  are  lift  (FI),  drag  (Fd), 
gravity  (Fg),  and  the  cross-stream  component  of  gravity  (Fgx).The 
gravitational  force  is  proportional  to  the  cube  of  the  particle  diameter, 
and  the  fluid  forces  are  proportional  to  the  square  of  the  diameter.  The 
result  is  a  sorting  mechanism.  (Adapted  from  Dietrich,  1988) 


Authors  using  the  former  method  can  control  the  secondary  currents 
and  the  bed  topography  using  these  equations,  but  one  needs  to  ask  whether 
the  simplifying  assumptions  made  in  order  to  solve  the  governing 
equations  are  justifiable  for  flow  in  bends. 

The  first  important  contribution  was  from  a  Dutch  engineer,  L  Van 
Bendegom,  in  1947,  who  formulated  a  model  similar  in  principle  to  the 
fundamental  model  described  in  the  previous  section,  and  was  the  first  to 
hypothesize  that  transverse  bed  slope  bears  a  simple  inverse  relation  to 
grain  size.  A  complete  appraisal  of  his  work  can  be  found  in  Allen  (1978). 
Independent  work  by  Rozovskii  (1957)  produced  very  similar  results. 
Rozovskii  developed  c  model  for  two-dimensional  flow  using  the  equations 
of  motion.  He  later  extended  his  analysis  to  consider  flow  in  three 
dimensions  and  finally  used  his  model  to  consider  the  development  and 
decay  of  the  helical  flow  cell.  Much  later,  Engelund  (1974)  used  an 
analysis  of  fluid  motion  developed  by  Rozovskii  to  produce  a  model  for 
flow  and  bed  topography.  He  claimed  that  the  velocity  defect  law  best 
describes  velocity  distribution  in  rivers.  He  first  approximated 
two-dimensional  bend  flow,  and  subsequently  included  some  second-order 
calculations  to  take  into  account  the  effect  on  the  flow  field  of  radial 
variations  of  depth  and  velocity.  Bridge  (1978)  used  Engelund's  approach 
for  bends  approximating  sine-generated  paths.  In  his  model,  there  is  no 
cross-stream  discharge  of  sediment  at  equilibrium,  and  grain  size 
increases  with  increasing  shear  stress,  so  that  the  maximum  sediment 
transport  is  towards  the  centre  of  the  channel. 
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There  has  been  much  criticism  of  the  approach  taken  by  Engelund  and 
Bridge.  Dietrich  and  Smith  (1984)  point  out  that  there  are  several 
deficiencies  in  the  flow  and  the  force  balance  equations.  In  particular, 
they  argue  that  convective  accelerations  should  not  be  considered 
second-order  effects,  but  are  important  enough  to  be  considered 
zero-order,  ie.  of  primary  importance.  Although  this  hypothesis  may  be 
correct,  Dietrich  and  Smith  found  that  their  model  is  very  sensitive  to  the 
convective  acceleration  terms.  This  means  that  to  apply  their  model 
successfully,  measurements  of  the  elevation  of  the  water  surface  have  to 
be  accurate  to  within  a  fraction  of  a  millimetre  (Anthony,  1987).  They 
claim  their  results  to  be  supported  by  the  flume  data  of  Yen  (1970).  The 
implication  of  their  argument  is  that  the  core  of  maximum  sediment 
transport  crosses  the  channel  in  a  bend.  Because  of  convective 
accelerations,  the  equilibrium  bed  slope  is  one  that  causes  sufficient 
cross-stream  bedload  transport  against  the  inward  secondary  currents, 
that  the  zone  of  increasing  boundary  shear  stress  in  the  pool  is  balanced 
by  a  convergence  of  sediment  transport;  the  cross-stream  sediment 
discharge  is  not  zero.  In  addition,  they  describe  a  more  sophisticated 
sorting  procedure  for  sand-bed  channels  where  the  migrating  shear  stress 
field  causes  skewing  of  dune  orientation.  In  the  upstream  part  of  the  bend, 
where  the  maximum  shear  stress  is  near  the  inner  bank,  the  inner  bank  end 
of  dunes  will  migrate  forward  more  rapidly  than  the  outer  bank  ends,  so 
the  dunes  will  be  skewed  across  the  channel.  Troughwise  currents  in  the 
lee  of  the  dune  cause  inward  transport  of  fine  particles  while  larger  ones 
roll  outwards.  Near  the  bend  exit  the  shear  stress  peak,  having  crossed 
the  channel,  causes  the  outer-bank  end  of  the  dunes  to  migrate  more 
rapidly  than  the  inner  bank  end,  causing  troughwise  transport  towards  the 
pool  even  against  the  action  of  the  main  helix  (figure  27).  This  is  in 
addition  to  the  simple  force  balance  sorting  scheme  proposed  by  Wilson 
(1973)  which  still  operates  on  the  stoss  side  of  the  dunes.  They  argue 
that  only  where  the  shear  stress  field  does  not  vary  in  the  downstream 
direction  (at  the  downstream  end  of  a  bend)  will  a  particle  force  balance 


Figure  27  Sediment  pathlines  through  a  sand-bed  meander  bend.  (Adapted 
from  Dietrich  and  Smith,  1984a). 
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Smith  and  Mclean  (1984)  further  developed  Dietrich  and  Smith's 
(1983)  model,  and  successfully  tested  it  against  the  flume  data  of  Hooke 
(1975). 

Zimmerman  and  Kennedy  (1978)  proposed  a  model  for  transverse  bed  slope. 
They  considered  the  spatially  varying  centrifugal  force  as  exerting  a 
torque  force  on  the  water.  They  balanced  this  force  against  the  friction  of 
the  bed.  They  tested  their  results  against  flume  data  and  agreement  was 
generally  good  except  in  the  outer  bank  region.  However,  their  model  is 
only  applicable  for  fully-developed  flow  and  can  only  produce  a  linear 
transverse  bed  slope,  although  linear  transverse  bed  slopes  are  seldom 
found  in  river  bends. 

Another  important  distinction  can  be  made  between  models  that  deal 
with  fully-developed  flow  and  models  that  deal  with  developing  flow. 
Most  models  deal  with  fully-developed  flow,  but  Odgaard  (1986a  &  b) 

argued  that  an  understanding  of  the  role  of  developing  flow  is  critical  if 
erosion  and  depostion  in  river  bends  is  to  be  understood  fully.  He  proposed 
a  model,  based  on  a  solution  to  the  equations  for  conservation  of  mass, 
conservation  of  momentum,  and  lateral  stability  of  the  stream  bed,  that 
accounts  for  both  developing  flow,  and  convective  accelerations.  With 
flow  conditions  varying  in  the  downstream  direction,  the  governing 
equations  are  complex  and  difficult  to  solve,  so  models  in  this  class  need 
many  simplifying  assumptions.  Kalkwijk  and  DeVriend  (1980)  wrote  a 
model  for  two-dimensional  flow  in  order  to  simulate  flow  in  river  bends 
where  the  depth  is  small  compared  to  the  width,  and  the  width  is  small 
compared  to  the  radius  of  curvature.  They  found  good  agreement  between 
the  results  of  their  model  and  flume  data.  Discrepancies  were  found 
between  the  measured  and  theoretical  velocity  profiles  at  the  outer  bank, 
and  at  the  entrance  to  the  bend  where  the  theoretical  results  were  rather 
high.  DeVriend  and  Geldof  (1983)  developed  Kalkwijk  and  DeVriend’s 
(1980)  model  and  tested  it  using  data  from  two  bends  on  the  River 
Dommel,  Holland.  The  agreement  between  measured  and  computed  results 
was  generally  good,  but  the  model  did  not  work  well  for  the  bend  exits 
because  they  did  not  adequately  account  for  the  effect  of  secondary  flow. 

Chang(1984)  developed  a  model  (FLUVIAL-12)  applicable  to  curved 
alluvial  streams  with  non-erodible  banks,  and  able  to  simulate  stream  bed 
changes  during  a  given  flow.  FLUVIAL-12  is  unique  in  that  it  can  account 
for  changes  in  stage.  The  model  incorporates  the  major  effects  of  helical 
flow,  and  performed  well  when  tested  with  data  from  the  San  Lorenzo 
River. 

Models  Used 

Two  models  were  applied  to  predict  the  ratio  of  outer  bank  to  average 
velocities  in  bendways:  those  of  Bridge  (1982)  and  Odgaard  (1988).  These 
models  were  unusual  in  that  their  authors  were  willing  to  make  them  fully 
available  to  us  and  assist  us  in  their  application.  Most  modellers  do  not 
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release  their  models  in  this  way  for  a  variety  of  reasons.  Others  supplied 
copies  of  papers  reporting  their  models,  but  not  computer  codes.  The  task 
of  rewriting  entire  codes  for  extremely  complex  models  was  simply 
beyond  the  resources  available  to  this  project. 

The  details  of  the  models  may  be  found  in  papers  authored  by  Bridge 
(1982)  and  by  Odgaard  (1988).  Some  additional  parameters  were  required 
for  application  of  the  models.  These  are  listed  in  Appendix  B. 

Menu-driven  FORTRAN  programs  were  produced  by  a  research 
associate,  Dr  Andrew  Markham  on  the  basis  of  complete  codings  supplied 
by  the  authors.  A  disk  with  the  menu-driven  programs  may  be  found  at  the 
end  of  this  report,  in  Appendix  C.  These  programs  were  used  to  produce 
estimates  of  the  maximum  depth  averaged  velocity  over  the  toe  of  the 
outer  bank  at  each  natural  and  laboratory  channel  bend. 


Results  of  Model  Applications 

The  results  are  listed  in  Table  11.  Bridge's  model  failed  for  bends 
with  very  low  radius  to  width  ratios  (Rc/w  <  1),  which  caused  the  model 
to  crash.  Odgaard's  model  failed  in  many  more  cases.  The  problem  was 
that  in  long  bends  the  model  predicted  negative  water  depths  at  the  inner 
bank,  leading  to  its  crashing.  Further  work  on  the  model  has  so  far  failed 
to  resolve  this  problem,  although  a  solution  must  be  possible.  Probably,  it 
will  be  essential  to  work  directly  with  Prof.  Odgaard  to  solve  it. 

Analysis  of  Results 


The  results  are  plotted  as  observed  versus  predicted  outer  bank 
velocities  in  Fig.  28.  The  agreement  is  generally  quite  good,  although 
systematic  errors  are  apparent  in  both  graphs.  In  the  case  of  Bridge's 
model  these  can  be  largely  explained  by  a  plot  of  the  percent  error,  defined 
by: 


error  = 


(observed  V^-  predicted  V^} 


observed  V, 


x  100% 


(12) 


tee 


versus  the  Rc/w  ratio.  The  results  are  plotted  in  Fig  29.  Bridge's  model  is 
seen  to  give  generally  excellent  accuracy  for  bends  with  Rc/w  values 
greater  than  2,  but  to  be  unacceptable  for  tighter  bends  with  Rc/w  <  2. 
This  accords  with  the  theory  of  bend  flow,  the  results  of  the  data  based 
approach  and  Bridge's  own  guidelines  on  the  use  of  his  model.  More  data 
are  needed  for  Odgaard’s  model,  but  the  preliminary  results  show  that  it  is 
prone  to  underestimating  the  toe  velocity  by  up  to  about  30  percent.  This 
occurs  because  for  bends  with  coarse  sand  or  gravel  beds  the  model  does 
not  predict  the  bed  scouring  at  the  outer  bank  which  is  observed  in  nature. 
Outer  bank  depths  are  only  marginally  greater  than  centerline  depths,  and 
likewise  outer  bank  velocities  are  only  slightly  greater  than  the  centerline 
velocity.  However,  the  errors  do  not  increase  markedly  at  low  Rc/w 
values,  and  the  model  is  consistent  right  down  to  Rc/w  =  0.8. 
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Table  11.  Results  of  Model 


Observed  Vtoe 

Rc/w 

Bridge's  Vtoe 

t 

0.80 

2.87 

0.90 

2 

0.80 

0.82 

3 

0.60 

0.75 

4 

1.10 

0.88 

5 

1.35 

1.75 

1.70 

6 

0.80 

1.25 

0.45 

7 

0.74 

1.27 

1.25 

8 

0.70 

1.27 

0.92 

9 

0.69 

2.92 

0.70 

10 

1.05 

7.92 

0.95 

1 1 

0.95 

6.58 

0  99 

12 

1.03 

3.76 

1.02 

13 

0.80 

1.12 

1.06 

14 

0.93 

0.80 

15 

0.83 

7.83 

0.72 

16 

0.74 

1.94 

0.99 

17 

0.70 

4.70 

0.78 

18 

0.98 

3.80 

1.59 

19 

1.60 

4.84 

1.39 

20 

1.46 

17.90 

1.48 

21 

1.53 

16.23 

1.48 

22 

1.62 

10.00 

1.4? 

23 

1.55 

9.34 

1.56 

24 

1.69 

8.23 

1.70 

25 

1.81 

13.85 

1.62 

26 

1.55 

9.64 

1.63 

27 

1.55 

9.26 

1.54 

28 

1.67 

11.75 

1.66 

29 

1.54 

21.62 

1.58 

30 

1.69 

11.92 

1.84 

31 

0.55 

2.72 

0.57 

32 

0.61 

2.42 

0.50 

33 

0.75 

2.00 

0.67 

34 

1.60 

4.85 

1.77 

35 

0.52 

4.5 

0.50 

36 

0.55 

4.5 

0.56 

37 

0.60 

4.5 

0.59 

38 

0.46 

8.0 

0.46 

39 

0.48 

8.0 

0.48 

40 

0.35 

3.0 

0.29 

41 

0.35 

3.0 

0.29 

42 

0.55 

0.55 

43 

0.73 

0.70 

44 

1.36 

1.06 

45 

1.39 

1.05 

46 

0.73 

0.76 

47 

0.68 

0.76 

48 

0.42 

0.35 

Tests 


Error  (%) 

Odgaard  Vtoe 

Error  (%) 

-12.500 

0.52 

35.000 

0.56 

30.000 

0.42 

30.000 

-25.900 

0.67 

39.100 

43.800 

0.57 

28.700 

•68.900 

-31.400 

-1.400 

9.500 

-4.200 

0.971 

-32.500 

13.300 
-33.800 
-1 1.400 

0.66 

10.800 

-62.200 

0.93 

5.100 

13.100 

-1.400 

3.300 

9.300 
♦0.645 
-0.592 
10.500 
-5.200 

0.645 

0.599 

-2.600 

-8.900 

1.35 

15.600 

-3.600 

0.45 

18.200 

18.000 

0.45 

26.200 

10.700 

-10.600 

0.60 

20.000 

3.800 

0.40 

23.100 

-1.800 

0.47 

14.500 

1.700 

0.47 

21.700 

0.000 

0.41 

10.900 

0.000 

0.43 

10.400 

17.100 

0.30 

14.300 

17.100 

0.30 

14.300 

0.000 

0.49 

10.900 

4.100 

0.66 

9.600 

22.100 

1.06 

22.100 

24.500 

1.08 

22.300 

-4.100 

0.57 

21.900 

-11.800 

0.57 

16.200 

16.700 

0.37 

1 1.900 
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Fig.  28  Results  of  Model  Tests 


Analytical  Approach  -  Bridge's  Model 
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Error  (7.)  Error  (7.) 


Fig.  29  Error  Distributions  as  a  Function  of  Rc/w 


Analytical  Approach  -  Bridge’s  Model 
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Prediction  of  Outer  Bank  Velocity  in  Natural  Channel  Bendways 

On  the  basis  of  this  test,  Bridge's  model  can  be  relied  upon  to  predict 
the  outer  bank  velocity  in  a  bend  to  within  approximately  +/-15  percent 
for  bends  with  Rc/w  greater  than  2.  This  would  appear  to  make  it  a  strong 
candidate  for  adoption  as  a  design  method.  For  engineering  design 
purposes  it  might  be  desirable  to  introduce  a  factor  of  safety  of  1.15  to 
the  predicted  velocity  to  ensure  that  any  error  is  on  the  safe  side. 
Bridge's  model  should  definately  not  be  used  for  bends  with  Rc/w  <  2, 
where  it  is  prone  to  large  errors. 

Bridge's  model  is  available  in  menu  driven  format  for  IBM  PC 
computers  or  compatible  machines,  and  it  is  relatively  easy  to  use 
(Appendix  A). 

Odgaard's  model  could  not  be  fully  tested  as  it  crashed  for  many  of 
the  bends  studied.  For  the  bends  at  which  it  did  work,  errors  were 
between  5  and  40  percent,  for  bends  with  Rc/w  ratios  between  0.8  and  5. 
Further  work  is  needed  before  Odgaard's  model  could  be  recommended  for 
use  as  a  design  method,  but  the  fact  that  it  does  not  fail  for  short,  tight 
radius  bends  is  very  encouraging. 

Conclusions 

The  prediction  of  outer  bank  velocity  in  bendways  is  a  complicated 
problem.  This  velocity  is  affected  by  a  number  of  factors,  which  are 
themselves  closely  inter-related.  A  data  based  approach  seeks  to  develop 
a  predictive  equation  by  examining  the  relationship  between  independent 
and  dependent  variables  illustrated  by  observations  in  real  flows.  An 
analytical  approach  attempts  to  use  basic  physics  and  the  laws  of  motion 
to  derive  equations  describing  bend  flow.  It  must  be  remembered  though 
that  even  the  "analytical"  approach  involves  assumptions  concerning  the 
three  dimensional  nature  of  the  flow,  and  calibration  using  empirical  data. 

To  a  scientist  the  analytical  approach  is  preferable  philosophically, 
provided  that  the  underlying  theory  is  sound.  It  is  not  yet  agreed  that 
this  is  so  for  the  models  tested  here,  and  so  they  merit  use  only  if  they 
can  produce  acceptable  results.  This  appears  to  be  the  case  for  the  Bridge 
model. 


Recommendations 

It  would  be  foolish  to  base  a  new  design  method  on  the  results  of  a 
single,  relatively  small,  study.  Consequently,  it  is  recommended  that: 

1.  Follow-up  work  be  undertaken  to  verify  the  findings  reported  here.  In 
particular,  further  experimental  work  is  needed  in  laboratory 
channels  with  a  mobile  bed  and  static  banks  to  investigate  how  the 
pattern  of  bed  scour  affects  the  distribution  of  velocity  adjacent  to 
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2. 


the  bank. 

Field  data  on  shear  stress  distributions  on  the  outer  bank  nre 
essential  if  that  approach  to  riprap  sizing  is  to  be  sustained. 

3.  A  major  initiative  on  flow  in  very  tight  bends  of  Rc/W  <  2  is  needed. 
It  is  clear  that  results  from  longer  radiio  bends  cannot  be  back 
extrapolated  into  this  zone.  The  flow  in  ^uch  bends  can  impinge  on 
the  outer  bank  at  a  very  acute  angle.  This  links  the  problem  to  that  of 
impinging  flow  in  braided  channels. 

4.  A  better  definition  of  outer  bank  velocity  is  needed.  This  study  has 
shown  that  in  many  cases  the  maximum  flow  attack  on  the  outer  bank 
occurs  on  the  bank  face  some  distance  above  the  toe,  rather  than  at 
the  toe  itself.  This  is  not  accounted  for  in  the  present  design 
approach. 

5.  For  natural  channels  it  appears  that  different  equations  may  be 
appropriate  for  bends  with  straight  and  meandering  approach 
channels.  If  this  is  true  in  general,  it  has  important  implications  lor 
river  training  and  stabilization. 

6.  Several  of  the  equations  produced  here  deserve  further  testing  and 
consideration  as  design  methods. 

7.  On  the  basis  of  this  study,  the  model  developed  by  John  Bridge  of  the 
State  University  of  New  York  at  Binghamton  appears  to  predict  outer 
bank  velocity  to  within  +/-15  percent  tor  bends  with  Rc/w  >  2.  It  is 
therefore  recommended  as  a  possible  design  approach  to  riprap  sizirg. 
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SUMMARY 


Bend  scour  due  to  bed  erosion  by  curved  flow  leads  to  deep  pools  close  to  the  outer  bank  in 
meander  bends. 

The  strength  of  helical  flow  which  drives  bend  scour  is  known  to  be  a  function  of  the  geometry  of 
the  bend  and  the  hydiaulics  of  flow  through  the  bend.  Helical  flow,  fast  near  bank  velocities  and 
deep  scour  all  tend  to  de-stabilise  the  outer  bank  at  a  bend.  Bendways  are  often  the  sites  of  severe 
bank  attack  and  retreat  as  a  result.  Consequently,  a  great  deal  of  time  and  effort  is  spent  on  the 
stabilisation  of  outer  bank  lines  using  dumped  stone  and  riprap. 

For  a  bank  stabilisation  scheme  to  be  successful  it  must  guard  against  both  the  erosive 
velocity  (and  shear  stress)  of  the  flow  and  the  deep  scour  of  the  bed  adjacent  to  the  bank. 
Experience  shows  that  failure  is  just  as  likely  to  occur  due  to  under-estimation  and  allowance  for 
toe  scour  than  to  under-design  with  respect  to  velocity. 

In  an  earlier  report  the  authors  examined  the  methods  available  to  predict  near  bank  velocities 
and  bank  shear  stresses  at  the  outer  bank  in  bendways  (Thome  and  Abt,  1990).  In  this  report  the 
authors  go  on  to  consider  the  methods  available  to  predict  near  bank  scour  depth.  Both  analytical 
and  semi-empirical  approaches  are  considered.  The  analytical  approaches  use  existing  and  available 
bend  flow  models  developed  by  John  Bridge  and  by  Jacob  Odgaard.  The  empirical  approach  is 
based  on  a  statistical  analysis  of  hydrometric  data  from  the  Red  River  in  Arkansas  and  Louisiana 
undertaken  by  the  first  author  on  behalf  of  the  Vicksburg  District,  US  Army  Corps  of  Engineers. 

These  predictive  techniques  are  tested  using  a  data  set  for  over  250  bends  assembled  from  a 
variety  of  rivers  around  the  world.  All  flows  referred  to  in  the  set  are  high,  in-bank  flows 
corresponding  to  ’formative’  discharges  in  the  channel.  The  data  set  covers  a  wide  range  of  sizes 
of  river  from  laboratory  channels  to  very  large  alluvial  streams.  It  also  encompasses  rivers  with 
bed  materials  ranging  from  sand  to  boulders.  The  rivers  display  both  freely  migrating  bends  with 
outer  banks  formed  in  easily  erodible  alluvium,  constrained  bends  with  resistant  outer  banks  and 
stable  bends  with  revetments  at  the  outer  bank. 

The  results  confirm  that  there  is  a  close  relationship  between  bend  geometry  and  scour  depth 
that  may  be  characterised  by  a  dimensionless  graph  with  the  ratio  of  bend  radius  of  curvature 
divided  by  crossing  width  (Rc/w)  on  the  x-axis  and  the  ratio  of  bend  maximum  scour  depth  to  mean 
crossing  depth  (dmax/dbar)  on  the  y-axis. 
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Application  of  the  bend  flow  models  produced  mixed  results.  The  Bridge  model  produced 
wide  scatter  and  generally  tended  to  over-estimate  scour  depth.  The  model  worked  better  for  long 
radius  bends  and  errors  increased  alarmingly  as  the  Rc/w  value  decreased  to  about  2.  It  is 
recommended  that  application  of  this  model  be  restricted  to  bends  with  Rc/w  greater  than  4  and  that 
it  be  recognised  that  errors  of  +50%  are  common  and  +100%  are  possible.  Bridge  has  developed  a 
new  model  which  should  do  better  and  will  supply  this  for  public  use  shortly  (Bridge,  personal 
communication,  1992). 

Odgaard’s  model  consistently  under-estimated  scour  depth.  For  smaller  rivers  predicted 
values  were  about  -50%,  but  on  larger  rivers  the  model  crashed.  However,  unlike  Bridge's  model, 
the  results  were  consistent  even  for  the  tightest  bends  and  actually  did  quite  well  for  Rc/w  =  1. 
Under-predictions  seemed  to  be  due  to  the  difficulty  the  model  had  in  correctly  predicting  the 
entrainment  and  transport  of  any  sediments  coarser  than  fine  sand.  If  this  problem  could  be  over¬ 
come,  the  model  has  great  potential.  It  should  be  noted  that  the  version  of  Odgaard's  model  used 
here  is  being  replaced  by  an  up-dated  and  improved  model  (Odgaard,  personal  communication, 
1992). 

Overall,  the  empirical  method  produced  the  best  agreement  between  observed  and  predicted 
scour  depths.  Practically  all  predictions  fell  within  +/-  50%  of  observed  values  across  the  whole 
range  of  river  scales,  scour  depths,  bed  materials  and  bend  geometries.  The  great  majority  of  the 
predictions  fall  within  a  band  of  +30%  to  -25%,  which  is  close  to  being  acceptable  for  engineering 
purposes. 

It  is  therefore  recommended  that  the  empirical  equation  be  further  tested  and  evaluated  for  use 
as  an  aid  to  scour  prediction  in  bends  where  more  sophisticated  methods  and  models  are 
unavailable  or  inapplicable.  This  equation  is: 

(dmax/dbar)  =  2.07  -  0.19  In  (Rc/w) 

where, 

dmax  -  maximum  scour  depth  in  bend 
dbar  -  mean  depth  at  the  upstream  crossing 
Rc  =  bend  radius  of  curvature 
w  =  width  at  upstream  crossing 

The  equation  refers  to  maximum  scour  during  high,  steady,  in-bank  flows.  It  cannot  be  used 
to  predict  low-flow  scour,  or  scour  variation  with  changing  stage  during  a  hydrograph.  Also,  care 
must  be  taken  if  the  equation  is  used  to  predict  scour  associated  with  flows  greater  than  bankfull 
stage.  If  there  is  significant  inter-action  between  the  in-channel  and  overbank  portions  of  the  flow 
then  scour  patterns  and  amounts  are  not  predictable  on  the  basis  of  in-bank  hydraulics  alone. 
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MAIN  TEXT 


Introduction 

Most  alluvial  rivers  have  a  meandering  planform.  Such  rivers  naturally  migrate  back  and 
forth  across  their  flood  plain  by  a  combination  of  relatively  orderly  meander  loop  growth  and 
downstream  progression  which  is  interupted  occasionally  through  abrupt  by-passing  of  acute  bends 
by  chute  and  neck  cut-offs.  Generally,  meander  growth  and  progression  occur  through  retreat  of 
the  outer  bank  and  advance  of  the  inner  bank  of  a  bend,  although  in  particular  cases  this  pattern 
may  be  reversed.  The  resulting  flood  plain  deposits  vary  widely  in  their  composition  and 
engineering  properties,  depending  on  the  depositionaiy  environment  in  which  they  were  produced. 
Easily  eroded  meander  belt  alluvium  generated  by  point  bar  deposition  at  the  inner  bank  in  bends  is 
interspersed  by  tougher  back  swamp  deposits  from  slack  water  areas  and  resistant  clay  plugs  in 
abandoned  channels  (Fisk,  1943;  1947). 

A  great  deal  of  engineering  work  is  undertaken  each  year  on  rivers  of  all  scales  from  small 
creeks  up  to  the  Lower  Mississippi  to  curtail  and  control  the  lateral  activity  of  meanders.  This  is  an 
essential  component  of  projects  to  improve  navigation,  increase  flood  capacity,  stabilize  banklines, 
decrease  flood  plain  destruction  and  reduce  the  downstream  dredging  requirement  of  the  river. 
Bank  stabilization  is  achieved  using  structures  made  from  a  wide  variety  of  materials  and 
combinations  of  materials  and  treatments.  For  example,  on  the  Mississippi  River  upper  bank  riprap 
and  lower  bank  articulated  concrete  mattress  (ACM)  are  used  to  protect  against  the  erosive  attack  of 
the  near  bank  flow,  together  with  bank  regrading  and  sub-surface  drainage  control  to  prevent  mass 
failure.  On  smaller  rivers  riprap  may  be  used  alone,  or  increasingly,  its  use  in  combination  with 
low  cost  alternatives  and  vegetation  is  a  common  solution. 

Bank  stabilization  assumes  particular  significance  where  a  levee  or  an  important  piece  of 
infra-structure  is  set  back  only  a  short  distance  behind  the  bank  line.  At  such  locations  retreat  of  the 
bank  cannot  be  allowed  because  it  would  put  the  levee  or  structure  in  jeopardy.  Yet  despite 
awareness  of  the  potentially  damaging  impacts  of  failure,  revetments  do  still  occasionally  fail, 
resulting  in  retreat  of  the  bank  line,  destruction  of  areas  of  the  flood  plain  and,  in  some  cases,  the 
loss  of  a  section  of  levee.  Riprap  failure  due  to  under-design  of  the  size  or  thickness  of  the  riprap 
blanket  are  comparatively  rare.  More  often  the  cause  of  failure  is  scouring  of  the  bend  pool 
adjacent  to  the  outer  bank  to  a  greater  depth  than  that  allowed  for  in  the  design  of  the  revetment. 
Scour  below  the  toe  of  the  revetment  may  trigger  mass  failure  of  the  whole  bank  (including,  or 
followed  by,  the  levee  above)  by  one  a  variety  of  mechanisms.  The  mechanisms  identified  as  being 
the  most  critical  include:  rotational  slip,  slab-type  collapse  and  retrogressive  flow  failure  (Turnbull 
et  ai.,  1966;  Toney,  1988).  Consequently,  it  is  important  to  be  able  to  predict  the  likely  maximum 
scour  depth  in  a  bend  accurately  when  designing  revetment 
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However,  past  studies  have  demonstrated  the  difficulty  of  making  accurate  predictions  of  the 
likely  magnitude  and  location  of  the  maximum  scour  depth  in  a  bend.  For  example,  some  studies 
have  indicated  that  the  severity  of  outer  bank  pool  scour  in  revetted  bends  actually  increases 
compared  to  free,  alluvial  meanders  (Friedkin,  1945;  Thome,  1988).  Other  studies  have 
concluded  that  the  evidence  is  equivocal  and  have  gone  so  far  as  to  suggest  that  stabilizing  the  outer 
bank  with  riprap  may  lead  to  a  decrease  in  scour  depth  (Harvey  and  Sing,  1989).  It  is  therefore 
vital  that  improved  approaches  to  scour  depth  prediction  in  meander  bends  be  developed  so  that 
scour  potential  may  be  accurately  predicted  and  allowed  for  in  the  design  of  the  revetment. 
Predictions  should  and  perhaps  must  allow  for  and  include  any  increase  in  outer  bank  pool  scour 
that  might  follow  stabilization  of  the  bank  line. 

A  further  problem  in  scour  pool  prediction  lies  in  the  fact  that  pool  scour  is  discharge  related. 
As  stage  rises  it  is  generally  found  that  pools  in  bends  are  scoured  deeper,  with  the  eroded  material 
being  deposited  on  the  shoals  at  the  crossings  between  bends.  On  the  falling  stage  at  the  end  of  the 
event,  material  is  re-eroded  from  the  crossings  and  deposited  in  the  next  bend  pool  downstream. 
Hence,  there  is  no  clear  indication  of  the  high  stage  bed  topography  once  the  flood  has  passed. 
Observations  of  scour  depth  at  low  flows  are  relatively  straight  forward.  However,  it  is  much 
more  difficult  to  make  observations  of  scour  depths  at  high  flows  and  so  reliable  field  data  arc 
scarce.  Yet,  it  is  usually  during  high  flows  that  maximum  scour  occurs  and  it  is  during,  or 
immediately  after,  high  flows  that  most  failures  occur.  This  explains  why  the  design  of 
stabilization  works  is  based  on  a  flow  of  relatively  long  return  period  such  as  the  10  or  20  year 
flood.  In  this  regard  it  is  desirable  that  methods  be  established  to  extrapolate  observations  or 
calculations  of  scour  depth  for  low  or  intermediate  flows  up  to  higher  flow  levels. 


Objectives 

In  view  of  the  numerous  problems  of  bank  instability  encountered  and  the  great  expenditure 
of  funds  on  bank  stabilization  it  is  desirable  that  improved  methods  be  developed  to  predict  'he 
ultimate,  maximum  scour  depth  in  bends  of  the  meandering  rivers.  Ideally  these  aspects  of  channel 
geometry  should  be  predicted  analytically,  from  the  distribution  of  flow  velocity  and  boundary 
shear  stress  in  the  bend  in  general,  and  close  to  the  outer  bank  in  particular.  But  at  present  our 
understating  of  flow  processes  in  bends  is  incomplete  and  our  ability  to  simulate  the  flow  and 
sedimentary  processes  numerically  is  limited.  Also,  existing  approaches  to  flow  and  sediment 
modeling  are  research  based  and  often  have  data  requirements  that  cannot  be  met  in  routine,  day-to- 
day  applications.  This  has  led  to  the  development  and  use  of  empirical  equations  and  curves  based 
on  field  and  laboratory  observation,  but  lacking  a  basis  in  the  theory  of  fluid  flow  and 
sedimentation  in  curved  channels. 
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The  technical  objective  of  this  project  is  to  examine  empirical  and  analytical  approaches  to 

scour  depth  prediction  in  the  light  of  increased  data  availability  through  river  and  flume  studies  and 

in  view  of  recent  advances  in  bend  flow  modelling.  The  specific  objectives  are: 

1 .  To  compile  an  extensive  and  reliable  data  base  on  bend  scour  depths  in  relation  to  bend 
geometry,  boundary  materials  and  flow  hydraulics; 

2 .  To  examine  the  applicability  and  accuracy  of  present  empirically  derived  predictors  of 
scour  pool  depth  and  location; 

3 .  To  test  the  applicability  and  accuracy  of  present  analytical  models  of  bend  flow  as 
predictors  of  scour  pool  depth  and  location; 

4 .  To  identify  if  the  erosion  resistance  and  mass  stability  of  the  eroding  outer  bank  plays 
any  identifiable  role  in  affecting  scour  pool  depth; 

5 .  To  establish  the  degree  to  which  scour  pool  depth  in  revetted  bends  increases  relative  to 
similar,  unrevetted  bends; 

6.  To  examine  the  difficulties  of  extrapolating  predictions  based  on  flows  of  short  return 
period  to  long  period  or  'design  flow’  conditions;  and 

7 .  To  recommend  a  best  approach  to  scour  pool  depth  and  location  prediction  for  alluvial 
and  revetted  bends  of  natural  channels  for  future  research  and  development. 


Review  of  Possible  Approaches  to  be  Adopted 
Bend  Migration  and  Bank  Failure 

Bend  migration  occurs  through  bank  retreat  and  usually  takes  place  as  a  result  of  erosion  of 
the  bank  and  bed  adjacent  to  the  bank  by  the  flow,  coupled  with  periodic  mass  failures  of  the  bank 
unde.-  gravity.  Processes  and  mechanisms  of  bank  erosion  and  retreat  have  been  reviewed  in  detail 
elsewhere  and  are  not  exhaustively  reported  here  (Thome,  1978;  US  Army,  1981;  Thome,  1982; 
Thome  and  Abt,  1989;  Thome  and  Abt  1990). 

During  failure  the  slump  debris  falls,  slides,  or  flows  to  the  toe  of  the  bank  under  gravity.  It 
is  removed  from  there  by  the  river  over  a  period  of  days,  weeks  or  months.  This  is  the  basal  clean¬ 
out  phase  of  bank  erosion.  While  in  place  at  the  toe  of  the  bank,  slump  debris  tends  to  stabilize  the 
bank  by  protecting  the  intact  bank  from  further  flow  erosion  and  by  buttressing  it  against  mass 
failure.  Hence,  it  is  the  rapidity  with  which  it  is  removed  and  the  bank  re-eroded  which  primarily 
determines  the  long-tenn  rate  of  retreat  of  the  bank  (Thome,  1978;  Lapointe  and  Carson,  1987). 
This  is  the  case  even  in  situations  where  the  most  obvious  mechanism  of  bank  retreat  is  mass 
failure  rather  than  direct  entrainment  by  the  flow.  Removal  of  slump  debris  and  continued  bank 
erosion  depend  on  the  near  bank  sediment  balance  between  in-coming  sediment  from  uostream  and 
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out-going  sediment  downstream,  which  in  turn  depends  on  the  distributions  of  near  bank  velocity 
and  boundary  shear  stress  around  the  bend.  If  supply  and  removal  rates  are  balanced  then  the 
elevation  of  the  toe  is  constant  and  the  bank  retreats  by  parallel  retreat.  However,  if  the  scour 
potential  of  the  flow  near  the  bank  exceeds  the  supply  of  material  from  bank  erosion  and  failures, 
the  bed  is  scoured  to  make  up  the  deficit  and  elevation  of  the  bed  is  lowered.  Bed  degradation  has 
the  effect  of  further  destabilizing  natural  banks  and  thereby  increasing  the  volume  and  frequency  of 
mass  failures.  In  this  way  the  sediment  flux  at  the  toe  is  balanced  and  further  scouring  is 
suppressed.  But,  if  the  outer  bank  is  stabilized  by  a  structure,  mass  failures  do  not  occur  and  the 
bank  is  non-erodible.  In  this  case  the  ultimate  scour  depth  is  limited  either  by  the  resistance  to 
erosion  of  the  substrate  or  by  the  decrease  in  bed  shear  that  occurs  when  the  pool  approaches  the 
hydraulically  determined  maximum  for  the  given  bend  geometry.  If  no  resistant  materials  are 
encountered  and  the  bed  material  is  not  capable  of  armoring  by  selective  entrainment,  then  very 
deep  scouring  may  occur  before  the  bed  topography  fully  adjusts  to  the  imposed  shear  stresses  at 
the  bed. 

In  this  respect,  models  of  flow  and  sediment  processes  in  bends  have  highlighted  the  close 
links  between  bed  scour  and  bank  erosion  (Thome,  1978;  Lapointe  and  Carson,  1987;  Thome  and 
Osman,  1988;  Thome,  1991).  Investigations  show  the  locus  of  bank  erosion  and  retreat  to  closely 
follow  that  of  the  scour  pool  in  a  bend.  Partly  this  is  the  case  because  high  flow  velocities  attack 
both  the  bed  and  bank  simultaneously,  but  also  the  very  existence  of  a  deep  scour  hole  close  to  the 
bank  is  itself  a  destabilizing  factor  with  regard  to  mass  failure.  A  deep  scour  hole  increases  both 
the  bank  height  and  its  steepness,  both  of  which  reduce  the  factor  of  safety  with  respect  to  mass 
failure.  There  is  evidence  that  the  critical  height  for  mass  instability  may  represent  a  limiting  factor 
on  scour  depth  in  bends  with  banks  formed  in  weak,  alluvial  materials  (Thome  and  Osman,  1988). 
Conversely,  deep  scour  holes  are  often  associated  with  strong,  cohesive  banks.  This  is  particularly 
the  case  for  resistant  outcrops  of  backswamp  or  clay-plug  materials  in  otherwise  alluvial  bends. 
These  outcrops  often  have  significant  scour  holes  associated  with  them. 

These  concepts  and  arguments  explain  why  an  increase  in  scour  depth  may  follow 
stabilization  of  a  bank  through  the  construction  of  a  revetment.  Care  must  be  taken  to  properly 
allow  for  this  increased  scour  in  the  design  of  the  structure  in  order  to  guard  against  failure  by 
launching.  Conversely,  scour  should  not  be  over-predicted  as  this  leads  to  over-design  of  the 
structure,  poor  cost  effectiveness  and  wasted  money. 

The  depth  of  scour  is  also  discharge  dependent  At  low  flows  pools  at  bends  tend  to  fill  with 
material  washed  off  the  crossing  upstream.  At  medium  flows  pools  are  scoured  of  this  temporarily 
deposited  material,  while  at  the  highest  in-bank  discharges  the  bed  is  scoured  down  further  still. 
This  makes  it  difficult  to  predict  maximum  bed  scour  purely  on  the  basis  of  bed  topography 
established  from  low-flow  hydrographic  surveys.  Even  if  high  flow  soundings  are  available,  it  is 


11 


usually  still  necessary  to  extrapolate  because  the  data  refer  to  a  flow  lower  than  the  design  flood  for 
the  structure. 

Studies  of  bend  migration  and  bank  retreat  have  established  that  the  rate  and  distribution  of 
bank  erosion  may  be  related  to  two  fundamental  factors:  firstly  the  flow  velocity  adjacent  to  the 
bank  (for  example  Hasegawa,  Parker,  Ikeda;  see  Ikeda  and  Parker,  1991  for  a  state  of  art  review) 
and,  secondly,  the  depth  of  scour  adjacent  to  the  bank  (for  example,  Odgaard;  see  Odgaard,  1989 
for  a  recent  summary).  Many  researchers  see  these  as  alternative  approaches  and  argue  that  if  one 
of  the  factors  is  successful  then  the  other  must  be  wrong.  It  is  now  well  established  that  in  nature 
banks  may  retreat  either  because  of  flow  erosion  by  direct  entrainment,  or  because  of  toe  scour  and 
mass  collapse  followed  by  basal  clean-out.  Hence,  the  two  factors  of  near  bank  velocity  and  near 
bank  depth  are  complimentary  to  one  another  and  in  different  circumstances  either  one  may  be 
primarily  responsible  for  controlling  the  rate  and  distribution  of  bank  retreat.  In  this  respect  near 
bank  velocity  and  scour  depth  are  recognised  not  as  alternative  factors  in  controlling  bank  stability 
and  retreat,  but  are  complimentary  factors. 

Modelling  Bend  Processes 

To  translate  this  qualitative  understanding  of  bend  flow  and  sediment  processes  into  a 
quantitative  and  predictive  approach  depends  either  on  the  statistical  analysis  of  empirical  data 
defining  the  relevant  variables,  or  on  the  application  of  analytical  models  for  depth  averaged  flow  in 
bends.  The  empirical,  or  data-based  approach  must  be  based  on  experience  from  a  wide  range  of 
types  and  scales  of  river.  The  limits  to  the  data-base  set  limits  to  the  applicability  of  the  method  and 
to  the  reliability  of  the  predictions.  However,  while  it  is  desirable  to  use  as  large  a  data-set  as 
possible,  close  attention  must  also  be  paid  to  data  quality.  If  spurious  or  unreliable  data  are 
included,  these  gready  diminish  the  value  of  the  data-base.  Hence,  the  empirical  approach  is 
essentially  data-dependent  and  relies  on  wide  ranging  but  high  quality  data.  Conversely,  to 
understand  the  rationale  for  the  selection  of  the  particular  models  it  is  necessary  to  review  briefly 
the  basis  for  two-dimensional  bend  flow  modeling. 

The  complete  equations  of  motion  for  flow  in  bends  cannot  be  solved  analytically.  Provided 
that  certain  simplifying  assumptions  can  be  made,  then  numerical  solutions  are  possible  (Smith  and 
McLean,  1982).  However,  controversy  in  the  academic  world  centers  on  the  nature  of  the 
assumptions  which  are  tenable. 

The  question  the  validity  of  assumptions  is  best  addressed  through  consideration  of  the 
relative  importance  of  each  of  the  terms  in  the  equations  for  downstream  and  cross-stream  slopes  at 
meander  bends.  These  equations  have  been  derived  from  the  governing  equations  for  the  vertically 
averaged  motion  through  bends  by  Smith  and  McLean  (1984).  The  derivation  requires  that  certain 
assumptions  be  made,  the  main  one  being  that  shear  stresses  due  to  lateral  boundary  layers  can  be 
neglected.  Dietrich  and  Whiting  (1989)  give  the  slope  equations  so  derived  as: 
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where  S  =  downstream  water  surface  slope,  Sn  =  cross-stream  water  surface  slope,  (tzs)b  = 

downstream  component  of  total  boundary  shear  stress,  r  =  fluid  density,  g  =  acceleration  due  to 
gravity,  h  =  flow  depth,  <us  >  =  vertically  averaged  downstream  component  of  fluid  velocity,  N  = 

n/R  (n  =  cross-stream  coordinate,  following  the  channel  centerline  and  positive  towards  the  left 
bank,  R  =  radius  of  curvature  of  the  channel  centerline),  s  =  downstream  coordinate,  parallel  to  the 
channel  centerline  and  <un  >  =  cross- stream  component  of  fluid  velocity.  In  a  simplified  form  the 

equations  may  be  written: 

S  =  Si  +  S2  +  S3  (3) 

sn  =  snl  +  Sn2  +  sn3  +  sn4  (4) 

Expressed  in  words,  the  first  equation  states  that  the  downstream  slope  (S)  is  the  sum  of  the 
slope  components  due  to  the  downstream  bed  shear  stress  (Si),  the  downstream  change  in 
momentum  (S2)  and  the  cross-stream  change  in  momentum  (S3).  The  second  equation  states  that 
the  cross-stream  slope  (Sn)  is  the  sum  of  slope  components  due  to  the  cross-stream  shear  stress 
(Sni),  the  centrifugal  acceleration  (Sn2),  the  cross-stream  change  in  momentum  (Sn3)  and  a  small 
term  produced  by  substitution  of  the  continuity  equation  into  the  equations  for  fluid  motion  (Sn4). 
Since  Sn4  is  small  it  can  safely  be  neglected  (Dietrich  and  Whiting,  1989). 

These  full  equations  show  that  the  actual  downstream  slope  cannot  always  be  approximated 
by  the  Darcy- Weisbach  equation  which  only  yields  S 1  and  which  is  only  applicable  to  uniform 
flow  (where  S2  and  S3  are  zero),  and  that  the  actual  transverse  slope  at  a  bend  does  not  depend 
only  on  the  centrifugal  acceleration,  but  is  also  affected  by  cross-stream  shear  stress  and  changes  in 
cross-scream  momentum.  A  complete  analysis  of  these  factors  is  beyond  the  scope  of  this 
proposal,  but  a  few  points  may  be  made. 

With  regard  to  the  downstream  slope,  S,  it  is  accepted  that  both  Sj  and  S2  are  important  but 
two  schools  of  thought  exist  regarding  the  significance  of  S3  .  Members  of  the  first  school  base 
their  approach  on  the  seminal  analysis  by  Engelund  (1974)  and  hold  that  the  S3  term  is  negligible. 
Notable  papers  by  Bridge  (1977;  1984)  and  by  Odgaard  (1986;  Odgaard  and  Bergs,  1988)  follow 
this  approach  and  present  field  and  flume  data  to  support  their  case.  The  second  school  base  their 
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approach  on  a  re-examination  of  the  equations  of  fluid  motion  undertaken  at  the  University  of 
Washington  in  the  late  1970s  by  a  team  of  researchers  including  Dietrich,  Dunne,  McClean  and 
Smith  (see  Dietrich  et  al.,  1979;  Dietrich,  1982;  Dietrich  and  Smith,  1983;  Dietrich  et  al„  1984; 
Smith  and  McClean,  1984).  They  hold  that  the  S2  and  S3  terms  are  of  comparable  magnitude  in 

most  natural  rivers  which  have  point  bars  and  in  which  the  local  radius  of  curvature  changes 
continuously  through  the  bend,  although  the  S3  term  may  well  be  small  in  laboratory  flumes  which 

lack  point  bars  and  have  bends  of  constant  radius.  With  respect  to  contrasts  between  sand-bed  and 
gravel-bed  rivers,  in  a  recent  comparison  reported  by  Dietrich  and  Whiting  (1989)  it  was  concluded 
that  the  ratio  of  S3  to  S2  in  gravel-bed  streams  was  similar  to  values  in  sand-bed  streams  despite 

marked  differences  in  curvature  and  bed  grain  size.  This  was  taken  to  demonstrate  that  bed 
topography  effects,  and  particularly  convective  acceleration  terms  due  to  topographic  steering  by 
the  point  bar,  are  the  dominant  controls  on  downstream  slope  at  bends. 

With  regard  to  the  cross-stream  slope,  Sn  ,  Dietrich  and  Whiting  (1989)  used  detailed  field 
observations  to  show  that  the  cross-stream  momentum  transfer  term,  Sn3«  is  only  negligible  if  flow 
at  the  bend  entrance  is  skewed  so  that  the  maximum  velocity  is  already  near  the  outer  bank.  In 
most  rivers,  where  consecutive  bends  are  of  opposite  curvature,  this  will  not  be  the  case  and  the 
cross-stream  discharge  and  momentum  flux  associated  with  the  movement  of  the  core  of  maximum 
velocity  from  near  the  inner  to  near  the  outer  bank  will  be  an  important  factor  in  determining  the 
cross-stream  water  surface  slope  in  both  sand-bed  and  gravel-bed  rivers. 


Approaches  Adopted 

On  this  basis  it  would  appear  that  models  based  on  the  Engelund  approximations  do  neglect 
significant  terms  and  that  models  based  on  the  full  equations  of  Smith  and  McLean  are  preferable. 
However,  derivatives  of  the  Smith  and  McLean  analysis  place  such  heavy  requirements  on  the 
accuracy  and  availability  of  data  that  they  cannot  be  applied  to  field  situations.  For  example,  to 
properly  account  for  the  effects  of  convective  accelerations  in  bend  flow  it  is  necessary  to  know  the 
water  surface  topography  in  the  bend  to  millimeter  accuracy  (Anthony,  1987;  Dietrich  and  Whiting, 
1989).  This  would  be  possible  to  measure  at  normal  flow,  in  small  rivers,  although  only  at  very 
great  expense.  However,  to  use  such  models  to  predict  the  longterm  bed  configuration  in  large 
rivers  would  not  be  possible  and  to  deal  with  high  flows  of  around  bankfull  discharge  in  even 
medium  sized  watercourses  is  not  feasible  with  present  technology. 

Recently,  in  an  evaluation  of  several  bend  flow  models  Markham  (1990)  concluded  that 
models  based  on  the  Engelund  approach  and  neglecting  the  convective  acceleration  terms  can  give 
reasonable  approximations  of  the  main  flow  and  morphological  features  of  meander  bends.  He 
found  that  the  main  features  of  bends  are  predictable,  but  that  in  detail  every  bend  is  different.  This 
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is  the  case  because  each  bend  is  the  unique  product  of  the  sequence  of  flow  events,  approach 
channel  orientation  and  sedimentary  materials  encountered  in  the  bed  and  banks.  These  factors  are 
stochastic  and  cannot  be  predicted  deterministically.  However,  for  practical  and  engineering 
purposes  there  is  no  justification  for  attempting  to  predict  them,  since  the  fine  details  of  flow  and 
morphology  are  usually  of  no  consequence  to  engineering  schemes  to  stabilize  the  bend.  This 
conclusion  has  also  been  reached  by  even  the  strongest  advocates  of  the  most  complex  bend  flow 
models,  who  agree  that  while  "blue  skies"  and  academic  studies  should  strive  to  be  geophysically 
correct,  much  simpler  models  may  be  safely  applied  for  practical  and  design  work  (Dietrich  and 
Whiting,  1989). 

The  most  promising  approach  is,  therefore,  to  use  the  latest  versions  of  models  developed 
from  the  Engelund  analysis  and  having  data  requirements  which  can  realistically  be  met  in  real 
world  situations.  In  a  previous  study  of  outer  bank  toe  velocities  in  bends  Thome  and  Abt  (1990) 
applied  PC  compatible  versions  of  models  by  Bridge  and  by  Odgaard  which  fulfil  these  criteria. 
They  were  found  to  give  outer  bank  velocity  predictions  which  were  better  than  those  obtainable 
from  empirical  curves  provided  that  close  attention  was  paid  to  the  limits  of  applicability  of  the 
models.  The  complete  findings  of  that  study  are  not  repeated  here  since  the  work  was  undertaken 
for  WES  under  contract  number  DACW39-89-K-0015  and  is  reported  in  detail  in  the  relevant 
document  (Thome  and  Abt,  1990). 

The  technical  approach  in  this  project  applies  and  tests  the  Bridge  and  Odgaard  models  as 
predictors  of  scour  pool  depth.  In  parallel  an  empirical  analysis  of  the  data  base  will  be  undertaken 
to  develop  morphometric  relationships  between  independent  variables  describing  the  bend 
geometry  and  dependent  variables  defining  the  location  and  severity  of  scouring  in  the  bend  pool. 
A  project  along  these  lines  was  recently  successfully  completed  using  data  from  a  1981 
hydrographic  survey  of  the  Red  River,  dealing  with  the  reach  between  Index,  Arkansas  and 
Shreveport,  Louisiana  (Thome,  1988).  Data  from  70  bends  was  used  to  produce  a  regression 
equation  for  the  prediction  of  maximum  bend  scour  (dmM)  on  the  basis  of  the  mean  depth  of  the 
approach  channel  at  the  crossing  upstream  of  the  bend  (dbar)  and  the  bend  geometry  represented  by 
the  ratio  of  bend  radius  (Rc)divided  by  width  at  the  upstream  crossing  (w).  The  regression 
equation  for  the  empirical  method  is: 

(dmax/dbar)  =  2.07  -  0.19  In  (R^/w)  (5) 

The  correlation  coefficient  was  0.8  and  was  statistically  significant  at  the  0.01  level  of 
confidence.  The  coefficient  of  determination  was  0.64,  indicating  that  variation  in  (Rc/w)  was  able 
to  account  for  64%  of  the  variation  in  (dmax/dbar)-  The  lower  limit  to  the  applicability  of  the 
equation  is  an  Rc/w  value  of  2.  This  is  consistent  with  the  observation  that  the  monotonic  increase 
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in  scour  depth  with  decreasing  Rc/w  seems  to  stop  at  this  value.  A  complete  account  of  the 
methodology  is  available  in  Thome  (1988). 


Data-Base:  Sources  and  Features 

Sources  of  Data 

The  data  used  in  this  study  come  from  a  wide  variety  of  types  and  sizes  of  rivers,  located  in 
different  physiographic  regions  and  from  different  parts  of  the  world.  The  data  were  compiled 
from  existing  reports  and  papers  supplement  by  information  solicited  in  letters  sent  to  researchers 
known  to  be  actively  interested  in  bend  flow.  Different  studies  and  researchers  supplied 
substantially  different  types  of  data.  Hence,  before  embarking  on  the  detailed  analysis  of  the  data, 
it  is  necessary  to  review  each  of  the  data  sources  and  the  nature,  scope  and  features  of  the  material 
supplied.  This  may  seem  unnecessary  to  readers,  but  it  is  in  fact  crucial.  Before  using  any 
observed  data  to  derive  or  test  empirical  and  analytical  models  it  is  vital  to  identify  the  particular 
characteristics  of  the  data  that  may  substantially  impact  the  results. 

Data-Base  Characteristics 

Thome  and  Abt  (1990).  The  first  data  used  came  from  Thome  and  Abt  (1990)  and  were 
compiled  in  a  project  concerned  with  the  prediction  of  depth  averaged  velocity  over  the  toe  of  the 
outer  bank  in  meanders.  The  base  data  for  natural  riven  taken  from  Thome  and  Abt  (1990)  are 
shown  in  Table  1.  These  data  were  compiled  especially  to  allow  application  of  the  Bridge  and 
Odgaaid  models  and  so  all  of  the  required  parameters  were  already  recorded.  The  data  had  been 
screened  in  the  previous  study  and  found  to  be  of  high  quality  and  reliability.  But  the  previous 
application  to  the  prediction  of  near  bank  velocity  did  not  involve  knowing  the  maximum  scour 
depth  at  each  bend  and  so  no  record  of  the  observed  value  was  made  in  Thome  and  Abt’s  report. 
In  this  study  the  original  reports  were  re-examined  and  the  maximum  scour  depths  observed  at  each 
study  bend  were  extracted  and  added  to  the  data  base.  It  was  noted  that  as  the  surveys  were  based 
on  surveying  of  a  finite  number  of  cross-sections  (between  three  and  seven  per  bend)  it  could  not 
be  guarantied  that  the  recorded  maximum  scour  depth  was  actually  the  maximum  observed 
anywhere  in  the  bend.  Also  it  was  noted  that  measurements  corresponded  to  high,  in-bank  flows 
in  the  rivers  rather  than  to  'design  flows'  of  long  return  period.  The  additional  data  used  in  this 
study  are  also  listed  in  Table  1. 

The  consistency  of  the  data  was  checked  by  computing  dimensionless  parameters  of  bend 
geometry  and  scour  depth.  The  bend  geometry  is  represented  by  the  ratio  of  bend  radius  of 
curvature  (Rc)  to  width  at  the  crossing  upstream  of  the  bend  (w).  The  scour  depth  is  represented 
by  the  ratio  of  bend  maximum  scour  depth  (dmax)  to  average  depth  at  the  crossing  upstream  of  the 
bend  (dbar).  The  results  are  listed  together  with  other  calculated  parameters  in  Table  2  and  are 
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plotted  in  a  dimensionless  graph  in  Fig.  1.  This  graph  uses  (Rc/w)  on  the  x-axis  and  (dmax/dbar) 
on  the  y-axis.  Theories  of  bendflow  (Bagnold,  1960),  flume  studies,  engineering  experience  and 
observations  on  natural  rivers  have  shown  that  bend  scour  may  be  up  to  3  times  the  mean  depth  in 
the  approach  channel  and  that  the  peak  scour  usually  occurs  in  bends  with  Rc/w  in  the  range  2  to  3. 
Examination  of  Fig.  1  indicates  that  these  data  are  consistent  with  these  rules  of  thumb  and  on  this 
basis  there  is  no  reason  to  doubt  their  validity. 

Red  River  Hydrographic  Surveys.  The  second  and  third  sources  of  data  in  this  study 
are  hydrographic  surveys  of  the  Red  River  between  Index,  Arkansas  and  Shreveport,  Louisiana. 
Two  surveys  were  used,  those  of  1981  and  1969.  The  Red  in  this  reach  is  highly  mobile,  with 
rapid  bank  erosion,  bend  migration  and  planform  evolution.  Consequently,  the  bed  topography, 
bend  geometry  and  planform  configuration  in  1981  bears  very  little  relation  to  that  in  1969.  Due  to 
the  changes  occuring  between  these  two  dates,  it  may  be  concluded  that  the  surveys  yield 
essentially  independent  data  sets  on  the  relationship  between  bend  geometry  and  bed  topography. 
The  base  data  for  1981  are  shown  in  Table  3,  which  uses  data  taken  from  Thome  (1988)  together 
with  additional  data  collected  in  this  study.  The  depths  in  this  data  set  are  referenced  to  the  water 
surface  profile  for  the  two  year  flow,  which  was  identified  by  Biedenham  et  al.  (1987)  as  a  good 
guide  to  the  channel  forming  flow  for  this  river.  However,  no  measurements  of  the  discharge 
corresponding  to  the  two  year  flow  at  each  bend  were  recorded  by  Thome  (1988).  Information  on 
the  variations  in  the  volumetric  discharge  for  the  two  year  flow  along  the  study  reach  were  supplied 
by  the  Vicksburg  District,  US  Army  Corps  of  Engineers.  Also,  a  large  amount  of  additional 
information  on  the  meander  wavelength  and  channel  sinuosity  were  required.  The  first  PI, 
working  with  the  research  assistants  working  on  this  project,  used  maps  and  aerial  photographs 
supplied  by  the  Vicksburg  District  to  make  measurements  of  the  relevant  planform  parameters  for 
the  study  bends.  Calculated  parameters  used  in  this  project  are  listed  in  Table  4.  The  Red  River 
base  data  from  the  1969  hydrographic  survey  are  shown  in  Table  5.  with  the  derived  data  listed  in 
Table  6. 

The  dimensionless  scour  depth  versus  the  Rc/w  plots  for  the  Red  River  surveys  are  shown  in 
Figs.  2  and  3.  Both  graphs  show  the  expected  distribution,  with  relative  scour  depth  increasing  as 
Rc/w  decreases,  and  peaking  for  Rc/w  values  of  around  two.  While  the  maximum  relative  scour 
depths  for  Rc/w  =  2  to  3  are  substantially  higher  in  1969  than  those  for  1981,  the  overall  shape  of 
the  data  cloud  is  the  same.  It  may,  therefore,  be  concluded  that  the  Red  River  data  are  consistent 
with  existing  knowledge. 

British  Gravel-Bed  Rivers.  The  fourth  data  set  comes  from  a  study  of  stable  gravel-bed 
rivers  in  the  United  Kingdom  reported  by  Hey  and  Thome  (1986).  The  data  extracted  from  their 
paper  together  with  information  on  bend  geometry  obtained  in  this  study,  are  shown  in  Table  7. 
All  of  the  required  data  for  this  study  were  present  except  the  wavelength  and  radius  of  curvature 
for  each  individual  bend.  Wavelengths  and  bend  radii  could  not  be  measured  in  the  field,  as  were 
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the  other  parameters  in  this  Hey  and  Thome's  data  set,  but  were  estimated  from  available  1 :25,000 
topographic  maps.  The  scale  of  the  maps  and  size  of  some  of  the  channels  was  such  that  while 
measuring  wavelengths  was  straightforward,  estimating  the  bend  radius  was  difficult  on  the 
smallest  rivers.  Thus  the  accuracy  of  Rc  must  be  reduced  in  this  data  set 

It  should  also  be  noted  that  the  observed  maximum  bend  scour  depths  are  based  on  a  single 
cross-section  in  each  bend  and  may  not  necessarily  represent  the  actual  maximum  for  the  bend. 
Sections  were  located  in  the  field  specifically  with  the  intention  of  representing  the  maximum 
expression  of  the  pool  geometry,  however.  It  would  therefore  be  expected  that  true  maxima  be 
similar  to  the  observed  values,  but  they  could  be  a  little  larger  than  those  observed  in  some  cases. 

The  calculated  data  obtained  in  this  study  are  listed  in  Table  8  and  the  variation  of 
dimensionless  scour  depth  with  Rc/w  is  shown  in  Fig.  4.  This  shows  the  expected  distribution, 
with  dimensionless  scour  increasing  as  Rc/w  decreases  and  peaking  for  Rc/w  of  about  two.  The 
highest  values  of  relative  depth  are  rather  lower  than  those  in  the  other  data  sets,  consistent  with  the 
fact  that  single  sections  do  not  always  intersect  the  point  of  deepest  scour. 

An  alternative  explanation  for  the  lack  of  dmax/dbar  values  greater  than  about  2.5  relates  to 
the  geomorphological  and  sedimentary  features  of  these  rivers.  The  rivers  are  all  gravel-bedded 
and  display  a  coarse  surface  layer  on  the  bed.  Bed  material  is  known  to  be  transported  at  high  in¬ 
bank  flows,  but  only  at  relatively  low  rates.  The  coarse  bed  surface  is  called  a  mobile  armor. 
Many  engineers  and  scientists  belive  that  armoring  of  this  type  limits  the  depth  of  scour  because 
selective  winnowing  of  the  finer  material  leads  to  the  scour  hole  becoming  covered  by  a  lag  layer  of 
coarse  sediment.  In  such  cases  the  maximum  scour  depth  for  a  bend  would  be  less  than  that 
observed  in  a  similar  bend  with  a  sand  bed,  where  scour  potential  was  not  limited  by  armoring. 
Other  researchers  maintain  that  in  a  live-bed  with  an  adequate  supply  of  bed  load  from  upstream, 
scour  depends  more  on  flux  inbalance  between  sediment  supply  and  removal  and  that  particle  size 
does  not  significantly  influence  scour  depth  (Emmett  Laursen,  personal  communication  at  ASCE 
conferences). 

Geomorphologically,  these  British  rivers  are  in  a  post-glacial  environment.  The  bed  and 
flood  plain  materials  through  which  these  rivers  flow  were  laid  down  after  the  last  ice-age,  between 
about  1 1,000  and  6,000  years  ago.  At  that  time  melt  water  was  abundant  and  the  rivers  were  more 
powerful  than  they  are  today.  A  case  can  be  made  that  scour  potential  in  such  cases  is  limited  by 
the  fact  that  the  contemporary  river  has  inherited  its  channel  and  sediments  from  a  more  energetic 
predecessor.  This  inference  is  that  present  rivers  are  "underfu"  with  regard  to  the  geometric  and 
geomorphic  features  of  the  channel  and  that  they  cannot  adjust  the  channel  to  suit  the  current  flow 
regime.  However,  Hey  and  Thome  (1986)  report  regime  equations  for  the  British  rivers  which  arc 
very  similar  to  those  for  other  gravel-bed  rivers  known  to  be  in  adjustment  (from  the  USA, 
Canada,  New  Zealand)  and  so  this  seems  unlikely. 
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It  is  concluded  that  the  gravel-bed  river  data  are  reliable,  that  the  channels  are  adjusted, 
stable,  alluvial  features  and  that  any  contrasts  to  the  data  from  sand-bed  rivers  may  be  legitimately 
attributed  to  the  effects  of  the  coarse  bed  material  in  these  riven. 

River  data  from  other  Researchers.  The  fifth  data  set  consists  of  river  data  supplied 
by  various  other  researchers  in  response  to  letters  of  enquiry.  Letters  were  sent  to  66  individuals 
who  are  known  to  be  working  on  bend  flow.  Of  those  who  responded  positively,  only  four 
researchers  sent  data  that  were  sufficiently  complete  in  a  timely  fashion.  Several  data  sets  lacked 
information  vital  to  the  application  of  the  models.  Follow-up  letters  were  sent  but  replies  have  so 
far  not  been  received.  Other  data  sets  have  arrived  later,  and  too  late  for  inclusion  in  this  report. 
Many  researchers  did  not  respond  at  all  and  several  more  politely  declined  to  share  their  results  until 
they  had  exhausted  their  own  analyses.  The  overall  response  was  disappointing,  but  not 
unexpected. 

The  data,  listed  in  Table  9,  cover  a  wide  range  of  scales  of  flow.  At  the  largest  scale,  they 
include  twelve  bends  from  the  River  Ganges  in  Bangladesh.  Since  the  data  were  supplied  to  us  it  is 
not  possible  to  check  how  it  was  derived  in  great  detail,  hence  uncertainties  concerning  the  da^i  arc 
greater  than  for  the  other  sources.  However,  all  the  data  are  reported  to  refer  to  high,  in-bank 
flows  and  should  be  comparable  to  the  other  four  data-sets.  Calculated  data  derived  from  that 
supplied  are  listed  in  Table  10.  Plotting  the  dimensionless  scour  depth  versus  the  Rc/w  ratio  (Fig. 
S)  shows  the  expected  shape  to  the  data  cloud,  but  with  marked  differences  to  the  other  data-sets. 
First  it  is  noticeable  that  all  of  the  bends  have  rather  low  Rc/w  values.  The  longest  bend  has  an 
Rc/w  of  only  7  and  most  are  below  3.  In  view  of  this,  deep  bend  scour  is  to  be  expected.  Even 
so,  the  dimensionless  scour  depths  are  still  surprisingly  high.  In  one  case  the  maximum  bend 
scour  is  over  eight  times  the  mean  crossing  depth.  In  6  cases  it  is  more  than  four  times  the  mean 
depth.  These  figures  are  out  of  line  with  all  the  other  data  collected  in  this  study  and  this  must  cast 
doubt  on  their  validity.  Examination  reveals  that  the  twelve  highest  points  in  Fig.  5  represent  the 
twelve  bends  on  the  River  Ganges.  It  seems  likely  that  these  data  are  not  directly  comparable  to  the 
rest.  The  explanation  probably  lies  in  the  rather  low  mean  crossing  depths  of  the  Ganges  for  a  river 
of  such  great  discharge  (Q  =  29,500  cumecs).  Regime  equations  for  the  Indian  sub-continent 
would  suggest  a  crossing  width  of  about  3,000  metres  and  a  mean  depth  of  about  7  metres.  The 
observed  mean  values  for  the  twelve  bends  arc  7,238  metres  and  3.61  metres,  respectively.  From 
this  it  may  be  concluded  that  the  regime  of  the  lower  Ganges  is  not  consistent  with  that  of  a  single 
thread  meandering  channel.  The  channel  is  much  wider  and  more  shallow  and  is  in  fact  in  the 
transition  zone  between  meandering  and  braiding.  This  is  consistent  with  my  field  observations 
that  the  river  displays  elements  of  both  meandering  and  braided  planforms.  Hence,  the  data  are  not 
a  fair  test  of  methods  developed  for  single-thread  meandering  channels  in  dynamic  regime.  The 
data  points  corresponding  to  the  lower  Ganges  have  been  circled  to  identify  them  in  Fig.  5  and  are 
not  used  further  in  the  analysis. 
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Laboratory  Flume  Data.  The  sixth  and  final  data  set  comprises  data  from  laboratory 
flume  channels  with  mobile  bed  materials  that  have  been  used  to  simulate  flow  and  bend  scour  in 
river  ben  is.  These  data  are  included  for  interest  as  it  is  felt  that  real  river  data  present  a  better 
vehicle  than  flume  channels  for  tests  of  practical  scour  predictors.  However,  physical  modeling 
using  mobile  bed  sediments  demonstrates  that  flumes  can  simulate  the  morphology  of  real  rivers 
and  so  the  data  are  included.  The  data  are  listed  in  Tables  11  and  12,  and  are  plotted  as 
dimensionless  scour  depth  versus  Rc/w  in  Fig.  6. 

Conclusions 

The  data  set  prepared  in  this  study  contains  data  for  257  bends  on  natural  rivers  (excluding 
the  12  on  the  Lower  Ganges).  It  also  contains  8  points  for  laboratory  flumes  for  comparative 
purposes.  The  data  have  been  examined  carefully  and  (with  the  exception  of  data  for  the  multi¬ 
channel  Lower  Ganges)  are  believed  to  be  representative  of  high  in-bank  flows  in  rivers  with 
mobile  bed  materials  and  predominantly  self-formed  channels.  It  is  known  that  some  of  the  bends 
have  alluvial  banks  formed  in  erodible  and  unstable  sediments  while  others  have  banks  that  are 
formed  by  less  erodible  and  more  stable  backswamp,  clay  plug  and  lithified  materials.  Others 
bends  still  have  banks  which  have  been  artificially  stabilized  and  which  are  therefore  non-erodible 
and  fully  stable  with  respect  to  mass  failure.  This  information  will  be  used  later,  in  the  analysis  of 
the  role  of  outer  bank  materials  and  stability  in  affecting  scour  depth,  but  the  initial  analysis  is  based 
entirely  on  flow  mechanics. 


Analysis  of  Data 

Results  of  Scour  Depth  Predictions 

The  data  listed  in  Tables  1  to  9  were  used  to  apply  the  analytical  models  of  Bridge  and  of 
Odgaard,  and  the  empirical  model  of  bend  scour  developed  by  Thome  (1988)  for  the  1981 
hydrographic  survey  of  the  Red  River.  When  assessing  the  results  it  must  be  remembered  that  the 
empirical  model  must  fit  the  Red  River,  1981  data  reasonably  well,  since  it  was  developed  from 
those  data  The  other  5  data  sets  do,  however,  present  a  legitimate  test  of  the  empirical  equation. 

Observed  versus  predicted  plots  for  maximum  scour  depth  for  each  of  the  data  sets  are 
shown  in  Figs.  7  to  12.  There  is  considerable  scatter  evident  in  all  the  graphs.  Errors  in  predicted 
scour  depths  are  plotted  against  bend  geomeay  in  Figures  13  to  18. 

Examination  of  Results 

Thorne  and  Abt  (1990).  The  Bridge,  Odgaard  and  Empirical  methods  were  each 
applied  successfully  to  the  data  from  a  range  of  rivers  contained  in  the  Thome  and  Abt  data  set. 
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The  results,  in  Figures  7  and  13,  show  that  overall  the  empirical  method  did  best,  with  a  reasonable 
agreement  between  observed  and  predicted  scour  depths  and  errors  generally  less  than  +/-  30%  and 
always  less  than  +/-  50%.  Errors  tended  to  be  greater  for  tighter  bends,  that  is  at  lower  Rc/w 
values,  it  should  be  noted  that  the  tightest  bends,  with  Rc/w  less  than  2  cannot  be  analysed  using 
the  empirical  method. 

The  Bridge  model  produces  predictions  which  fall  around  the  line  of  perfect  agreement,  but 
which  involve  considerable  scatter.  Errors  are  most  pronounced  for  tighter  bends  and  rise  to 
greater  than  200%  for  the  most  acute  bends.  For  longer  bends  the  method  does  much  better,  with 
errors  generally  in  the  range  +/-  50%.  Poor  performance  in  short  radius  bends  was  a  feature  of  the 
Bridge  model  in  predicting  outer  bank  toe  velocity  (Thome  and  Abt,  1990)  and  so  this  result  is  not 
unexpected.  However,  even  for  the  longer  radius  bends  some  overestimates  of  80  to  90%  occur. 

The  Odgaard  model  underestimated  bend  scour  in  all  but  a  few  cases,  with  errors  as  great  as 
60%.  Examination  of  Figure  7  shows  that  predictions  are  quite  good  for  bands  with  scour  depths 
up  to  about  7  metres  deep,  but  that  they  seriously  under-estimate  the  observed  depths  for  deeper 
rivers.  This  problem  was  a  direct  result  of  the  difficulty  of  getting  the  model  to  run  for  large  rivers. 
Error  messages  and  warnings  associated  with  the  generation  of  negative  depths  at  the  inner  bank 
caused  problems  in  the  large  rivers.  Also,  the  model  did  not  deal  well  with  coarse  sediments.  If 
the  bed  material  was  coarser  than  sand  then  all  scouring  was  suppressed  and  the  predicted 
maximum  scour  depth  was  close  to  but  a  little  more  than  the  mean  depth  in  the  approach  channel. 
This  problem  was  only  overcome  by  using  a  sand  grain  size  even  when  the  indicated  bed  material 
was  coarser.  This  finding  corresponds  to  earlier  experience  with  this  version  of  the  Odgaard 
model. 

Red  River  Hydrographic  Survey  1981.  The  Red  River  survey  of  1981  provided  the 
data  for  development  of  the  empirical  (regression)  method  and  it  is  therefore  self-evident  that  the 
method  must  work  fairly  well  for  those  data.  Even  so,  the  predictions  are  quite  good,  with  errors 
generally  less  than  +/-  20%.  Errors  are  greater  for  the  tighter  bends. 

The  data  from  the  Red  River  include  many  tight  bends  which  led  to  scour  over-estimates  of 
up  to  100%  by  the  Bridge  model.  The  predictions  were  better  for  longer  radius  bends,  with  errors 
being  in  the  range  +  50%  to  -  20%  for  Rc/w  greater  than  4.  Unfortunately,  many  bends  in  practical 
problems  have  Rc/w  values  less  than  4. 

The  Odgaard  model  does  poorly  for  these  data,  all  which  come  from  bends  which  arc  really 
too  long  for  the  model  in  its  present  form.  The  model  systematically  under-estimates  the  observed 
scour  depth  with  errors  in  the  general  range  of  -40  to  -60%. 

Red  River  Hydrographic  Survey  1969.  The  second  set  of  data  from  the  Red  River 
produce  similar  results  to  the  first.  The  empirical  method  does  best,  with  predicted  values 
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straddling  the  line  of  perfect  agreement.  Errors  are  in  the  range  +/-  30%  over  quite  a  wide  range  cf 
R c/w  values. 

The  Bridge  model  is  again  prone  to  large  over-estimates  of  scour  depth,  especially  for  short 
radius  bends,  with  errors  uo  to  100%.  Performance  is  markedly  improved  for  less  acute  bends 
with  Rc/w  values  greater  than  4.  For  these  bends  the  Bridge  model  usually  under-estimates  the 
observed  scour  depth  by  about  20%. 

Odgaard's  model  under-estimates  scour  by  30  to  80%  across  the  whole  range  of  bend 
geometries.  This  is  consistent  with  its  performance  in  the  other  applications. 

British  Gravel-Bed  Rivers.  The  fourth  data  set  differs  from  the  others  in  dealing 
exclusively  with  gravel  and  cobble-bed  streams.  However,  despite  this  marked  contrast  in  bed 
materials,  the  empirical  method  continues  to  produce  the  best  agreement  between  observed  and 
predicted  scour  depths  (Figure  10).  Examination  of  the  error  distribution  (Figure  16)  shows  nearly 
all  errors  to  be  within  the  range  +50  to  -25%,  over  a  wide  range  of  bend  geometries. 

The  Bridge  model  did  not  work  well  for  the  gravel-bed  river  data.  It  over-predicted  scour 
depth  in  most  cases  (Figure  10),  with  errors  as  groat  as  300%  (Figure  16).  Similarly  to  the 
previous  data  sets,  the  most  serious  errors  occurred  in  the  tight,  short  radius  bends  with  Rc/w  less 
than  4  and  some  excellent  predictions  are  made  for  some  bends.  However,  even  for  longer  bends 
over-estimates  of  100%  still  occur  in  7  cases. 

The  Odgaard  model  was  not  applied  to  the  gravel-bed  river  data.  In  view  of  its  performance 
in  the  previous  tests,  it  was  clear  that  the  coarse  gravel,  cobble  and  boulder  sediments  common  to 
all  the  rivers  in  this  set  would  confound  the  model  and  result  in  predictions  of  minimal  scouring 
below  the  mean  depth  in  the  approach  channel. 

River  data  from  various  other  Researchers.  The  fifth  data  set  covered  a  variety  of 
rivers,  but  was  edited  to  remove  the  observations  from  the  Lower  Ganges  in  Bangladesh.  The 
examination  of  the  dimensionless  scour  depth  indicated  that  the  values  of  6  to  8  for  dmax/dbar 
differentiated  these  data  from  the  great  majority  of  the  other  observations  of  bend  scour  and  channel 
geometry  in  single-thread  meandering  rivers. 

The  plot  of  observed  versus  predicted  scour  depth  (Figure  1 1)  shows  that  the  points  for  the 
empirical  method  scatter  around  the  line  of  perfect  agreement  Predictions  are  better  for  the  deeper 
scours,  but  the  method  over-estimates  the  maximum  depth  in  some  of  the  shallower  bends,  by 
100%  in  one  case,  but  more  commonly  by  about  50%. 

The  Bridge  model  over-estimates  scour  depth  in  all  cases,  with  errors  of  at  least  50%  and 
greater  than  100%  in  four  cases. 
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The  Odgaard  model  does  quite  well  for  bends  with  scour  depths  of  less  than  6  metres,  but  it 
then  does  poorly  for  the  deeper  rivers.  This  performance  is  in  line  with  experience  in  the  previous 
applications. 

Laboratory  Flume  Data.  Application  to  the  laboratory  flume  data  showed  both  the 
empirical  and  Bridge  approaches  to  be  fairly  reliable.  Points  for  both  methods  scatter  around  the 
line  of  perfect  agreement  in  Figure  12,  with  errors  generally  in  the  range  of  +/-  20%  shown  in 
Figure  18.  However,  the  empirical  method  did  seriously  over-predict  scour  in  one  case  (by  60%) 
and  for  the  same  case  the  Bridge  model  also  did  comparatively  poorly,  being  off  by  33%.  This  is 
the  only  case  where  the  Bridge  model  actually  performs  slightly  better  than  the  empirical  method. 

Summary  data  for  Ail  Rivers 

The  separate  examination  and  plotting  of  the  data  from  the  different  sources  may  appear  long 
and  time-consuming,  but  it  was  essential  to  establish  the  consistency  of  the  results  across  the  range 
of  river  types,  researchers  and  measurement  techniques  employed  by  the  different  studies.  On  the 
basis  of  the  examination  of  the  results,  it  is  clear  that  the  performance  of  the  models  is  in  fact 
generally  consistent  across  all  the  data  sets  and  that  it  is,  therefore,  safe  to  compile  all  the  results 
together  in  order  to  produce  a  comprehensive  overview.  Summary  results  arc  shown  for  all  the 
data  in  Figures  19  and  20. 

Figure  19  compares  observed  and  predicted  scour  depths  for  all  256  bends.  The  data  cover 
maximum  scour  depths  that  range  from  a  few  centimetres  in  flume  channels  up  to  about  seventeen 
metres  in  large  rivers.  This  encompasses  the  likely  scour  in  all  but  the  world's  greatest  rivers. 
Wiih  regard  to  very  large  rivers,  data  for  the  Ganges  River  had,  unfortunately,  to  be  excluded 
because  of  their  inconsistency  with  the  other  observations.  Data  for  selected  bends  of  the  Lower 
Mississippi  have  already  been  examined  in  another  study  and  may  be  reviewed  there  (Thome  and 
Hubbard,  1991). 

Generally,  the  empirical  method  produces  the  best  overall  agreement.  The  points  cluster 
around  the  line  of  perfect  agreement  and  errors  appear  to  be  randomly  distributed.  Figure  20 
shows  that  practically  all  the  predictions  fall  with  +/-  50%  of  the  observed  values  across  the  whole 
range  of  scour  depths  and  bend  geometries,  and  the  vast  majority  fall  within  a  band  from  +30  to  - 
25%.  This  compares  very  favourably  to  the  analytical  predictions. 

Bridge's  model  produces  wide  scatter,  with  a  general  tendency  to  over-estimate  scour  depth. 
Examination  of  the  distribution  of  errors  a  function  of  bend  geometry  confirms  the  earlier  results  of 
Thome  and  Abt  (1990)  in  that  errors  increase  markedly  as  the  Rc/w  value  decreases  towards  a  value 
of  2.  It  has  been  recognised  that  in  such  tight  bends  the  flow  may  strike  the  outer  bank  at  a  steep 
angle,  driving  reversed  flow  upstream  of  the  apex.  In  this  respect,  attack  of  the  outer  bank  is 
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associated  with  impinging  flow,  and  this  has  found  to  be  an  intractable  problem  in  both  single 
thread  and  braided  rivers. 

Based  on  these  findings,  the  applicability  of  the  Bridge  model  is  limited  to  relatively  long 
radius  bends,  with  Rc/w  ratios  greater  than  about  4.  Even  in  these  bends  over-estimates  of  50%, 
and  occasionally  as  much  as  100%,  must  be  expected. 

Odgaard's  model  systematically  under-predicts  scour  depth.  For  the  smaller  rivers  points 
trend  around  a  line  of  -50%,  but  the  model  crashed  and  produced  little  further  increase  in  predicted 
scour  depth  once  the  observed  depth  passed  about  10  metres.  However,  Figure  20  shows  that 
Odgaard's  model  was  consistent  even  for  the  tightest  bends,  and  actually  did  well  for  cases  where 
Rc/w  was  around  unity.  If  the  problems  of  bed  material  mobility  and  application  to  large  rivers  can 
be  solved,  then  Odgaard's  model  has  real  potential,  especially  for  tight,  short  radius  bands  where 
the  empirical  model  is  inapplicable  and  the  Bridge  model  is  inaccurate. 

Conclusions 

This  study  has  assembled  a  fairly  large  and  internally  consistent  data-set  for  scour  at  bends. 
Application  of  three  different  approaches  to  scour  predicdon  has  illustrated  the  present  problems  of 
making  accurate  predictions  in  bends  of  different  geometry  and  size.  The  empirical  method 
produces  the  best  overall  agreement,  but  for  longer  radius  bends  Bridge;s  model  does  nearly  as 
well.  The  empirical  method  cannot  be  used  for  tight  bends  with  Rc/w  =  2  being  a  lower  boundary. 
Although  it  does  produce  predictions.  Bridge's  model  should  not  be  used  for  shorter  radius  bends, 
where  it  seriously  over-predicts  scour  depth.  Odgaard’s  model  remains  consistent  in  under¬ 
estimating  scour  depth  by  about  50%  over  a  wide  range  of  bend  geometries.  Accuracy  actually 
improves  for  very  tight  bends.  Hie  model,  therefore  shows  promise,  but  could  not  be  used 
routinely  in  its  1988  version. 

It  should  be  noted  that  both  Bridge  and  Odgaard  are  continuously  developing  and  upgrading 
their  models.  I  understand  from  both  researchers  that  the  problems  and  difficulties  noted  in  this 
study  are  consistent  with  their  experience  and  that  they  have  been  addressed  in  subsequent  editions 
of  the  models.  The  latest  versions  will  be  supplied  to  me  to  replace  the  earlier  ones  as  soon  as  they 
are  ready,  but  nothing  has  appeared  as  yet  (May  1992).  Updated  predictions  will  be  produced 
when  the  new  versions  arrive. 


Scour  in  Revetted  and  Free  Meanders 

It  was  discussed  in  the  Introduction  that  there  are  reasons  for  thinking  that  scour  depths  in 
revetted  bends  are  deeper  than  for  free  meanders  of  similar  geometry.  The  primary  purpose  of 
this  study  is  to  assess  the  capabilities  of  the  models  to  predict  scour  in  revetted  bends,  and  so  this 
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pan  of  the  analysis  concentrates  on  bends  with  revetted  outer  banks.  However,  as  the  data-set 
includes  both  revetted  and  free  meanders,  the  bends  were  split  along  those  lines,  to  allow 
comparison  of  revetted  and  free  conditions.  Bends  which  were  wholly  or  partially  constrained  by 
resistant  outcrops  of  clay,  backswamp  deposits  or  Pleistocene  materials  were  excluded.  Such 
bends  would  be  expected  to  have  scour  depths  intermediate  between  revetted  and  free  states,  and  so 
a  clearer  picture  should  emerge. 

The  observed  and  predicted  results  for  revetted  bends  are  shown  in  Figure  21.  They  are 
broadly  in  line  with  the  overall  study.  Bridge  s  model  produces  the  most  scatter  and  generally 
over-predicts  scour,  Odgaard's  model  systematically  under-predicts  and  the  empirical  model  does 
comparatively  well.  Examination  of  the  distribution  of  errors  (Figure  22)  shows  how  errors  in  the 
Bridge  model  increase  for  tight  bends.  His  model  looks  inapplicable  to  revetted  bends  with  Rc/w 
less  than  4.  But  for  longer  bends  the  model  does  quite  well,  with  predictions  in  the  approximate 
range  +15  to  -  25%.  The  empirical  model  is  successful,  but  as  a  considerable  number  of  the 
revetted  bends  come  from  the  Red  River,  1981  hydrographic  survey,  this  is  not  a  stringent  test  for 
the  method.  Errors  are  evenly  distributed  fro  a  range  of  Rc/w  values  and  fall  in  the  range  +25  to  - 
15%.  Odgaard's  model  seriously  under-predicts  scour  in  all  cases. 

In  the  Red  River  study  of  Thome  (1988),  a  separate  empirical  analysis  was  developed  for 
revetted  bends.  Like  the  general  equation  for  all  bends,  this  was  based  on  the  logarithm  of  ((Rc/w) 

-  2)  as  the  x  variable.  The  least  squares  regression  equation  was: 

(dmax/dbar)  =  2.15  -  0.27  In  ((Rc/w)  -  2)  (6) 

To  examine  if  this  equation  had  potential  for  the  expanded  data-set  a  semi-log  plot  of  (Rc/w) 

-  2)  versus  dmax/dbar  was  plotted  (Figure  23).  The  scatter  shows  a  linear  trend  to  the  data  with  a 
negative  slope.  Consequently,  the  data  are  consistent  with  the  form  of  the  empirical  revetment 
scour  equation.  Figure  24  shows  the  observed  and  predicted  scour  depths  for  the  empirical  and 
revetment  equations.  Figure  25  shows  the  distribution  of  errors  as  a  function  of  bend  geometry. 
There  is  no  obvious  improvement  in  prediction  for  the  revetment  equation,  and  it  actually  does  less 
well  in  the  revetted  long  radius  bends  of  the  Missouri  River.  On  the  basis  of  this  test  it  does  not 
appear  that  the  equation  for  revetted  bends  is  any  better  than  that  for  all  bends. 

Figure  26  shows  a  semi-log  plot  for  the  free  meanders.  Comparison  with  the  revetted  bends 
in  Figure  23  shows  that  the  two  data  clouds  overlap  to  such  a  degree  that  there  is  no  easy  way  to 
discriminate  between  them.  This  suggests  that  the  empirical  analysis  based  on  a  semi-log  equation 
is  equally  applicable  to  bends  of  all  outer  bank  types.  This  finding  is  in  contrast  to  the  findings  of 
the  earlier  study  on  the  Red  River.  There  are  at  least  two  explanations  for  this.  Firstly,  it  comes 
about  because  the  wider  range  of  conditions  in  the  rivers  of  the  large  data  set  have  swamped 
differences  due  to  bank  condition.  Thus,  while  bank  effects  may  be  identified  on  a  particular  river. 
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they  cannot  be  discerned  in  a  data-set  for  many  different  rivers.  Secondly,  the  parameters  used  in 
this  analysis  may  tend  to  collapse  data  for  different  bank  conditions  together.  This  is  the  case 
because  the  width  and  mean  depth  at  the  crossing  are  not  themselves  independent  of  the  bank 
condition.  Rivers  with  stiff  banks  tend  to  be  narrower  and  deeper  than  those  with  erodible  banks 
Hence,  in  using  crossing  width  and  mean  depth  to  scale  bend  geometry  and  scour  depth,  the  bank 
type  effects  are  being  implicitly  taken  into  account.  Consideration  of  the  impacts  of  revetments  on 
the  regime  geometry  would  help  to  resolve  this  issue,  but  is  getting  beyond  the  scope  of  this  study. 

It  can  be  concluded  though  that  the  empirical  method  constitutes  a  robust  predictor  of  scour 
depth  for  bends  of  a  variety  of  bank  types  and  conditions,  with  errors  usually  within  the  range  +/* 
30%. 
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Extrapolation  to  High  Flows 


The  flows  referenced  in  this  report  are  all  high,  in-bank  flows  that  are  close  to  'bankfull' 
discharge.  In  geomorphological  terms,  they  are  believed  to  be  the  channel  forming  flow  for  alluvial 
rivers:  that  is  the  flow  to  which  the  major  morphological  and  geometrical  features  are  adjusted. 
That  is  not  to  say  that  higher  flows  do  not  have  significant  impacts  on  the  channel  and  the  riparian 
zone.  It  is  well  known  that  scour  of  both  in-channel  pools  and  overbank  areas  can  be  considerable 
during  flood  flows.  However,  over  longer  periods  of  say  20  to  50  years,  the  impacts  of  high 
magnitude,  low  frequency  events  are  not  as  important  in  forming  the  landscape  as  the  impacts  of 
flows  at  and  around  bankfull  flow. 

In  the  case  of  revetment  design,  it  is  essential  that  the  structure  be  able  to  withstand  the  short 
term  impact  of  high  magnitude  events,  as  well  as  being  in  sympathy  with  the  hydraulic  geometry 
and  fluvial  morphology  of  the  river,  which  are  products  of  lesser  events.  Thus,  although  the  use  of 
flows  around  bankfull  have  a  role  in  determining  the  alignment  and  orientation  of  the  bank  line,  the 
strength  and  scour  protection  must  guard  against  short  duradon  but  higher  flows.  The  highest  flow 
to  be  guarded  against  is  termed  the  'design  discharge’  and  will  usually  be  a  flow  of  quite  long 
return  period.  Depending  on  the  level  of  acceptable  risk  and  the  value  of  the  land  being  protected, 
return  periods  of  20  or  even  50  years  may  be  used. 

In  teims  of  application  of  the  scour  predictors  tested  here  this  poses  a  problem  because  the 
flows  referenced  have  return  periods  of  the  order  of  1  to  5  years. 

In  order  to  solve  this  problem  it  is  necessary  to  know  something  about  the  flow  stage 
associated  with  the  design  flood  in  relation  to  the  elevation  of  flood  plain  and,  more  particularly,  the 
top  of  the  revetment.  In  general  terms,  the  stage  may  be  classed  as  either:  at  or  below  top  bank; 
slightly  overbank;  or  significantly  overbank. 

For  in-bank  or  bank  top  flows  little  difficulty  arises.  The  analytical  and  empirical  methods 
described  here  are  specifically  designed  to  deal  with  high  in-bank  flows  under  which  the  bed  is 
highly  mobile  and  to  which  the  geometry  of  the  channel  is  adjusted.  In  this  respect  the  design  flow 
should  be  ideal  for  application  of  the  methods.  Where  the  design  flood  does  not  overtop  the 
revetment  and  is  contained  within  the  channel,  no  extrapolation  should  be  necessary. 

Design  flows  which  are  slightly  out  of  bank  pose  more  of  a  problem,  since  by  strict 
definition  they  are  beyond  the  scope  of  the  models  and  methods  presented  here.  This  is  the  case 
because  once  the  flow  goes  out  of  bank  strong  eddies  develop  at  the  inter-face  between  the  channel 
and  flood  plain  elements.  These  eddies  induce  water  surface  topography,  momentum  exchanges 
and  shear  stress  distributions  which  cannot  be  predicted  on  the  basis  of  in-bank  flow  models  alone. 
However,  provided  that  the  depth  of  flow  on  the  flood  plain  is  small,  it  may  be  possible  to 
extrapolate  the  approach  for  in-bank  flows  up  to  out  of  bank  condith  -.is.  This  is  the  case  because 
relatively  high  flood  plain  roughness  and  low  flood  plain  flow  depth  combine  to  produce  low 


27 


velocities  of  flow  in  comparison  to  those  in  the  channel.  Hence,  the  overbank  area  may  be 
considered  more  as  a  reservoir  than  an  extended  channel.  Under  these  circumstances,  water  flow 
tends  to  follow  the  line  of  the  channel,  with  stable  boundaries  between  the  flowing  in-bank  and 
stadc,  overbank  water.  A  notional  tine  extending  up  from  the  bank  top  to  the  water  surface  may  be 
used  to  delineate  the  channel  for  modelling  purposes  and  the  methods  described  here  may  be 
applied  using  the  notional  geometry  of  the  cross-section  and  planform  of  the  flowing  channel. 
Experience  suggests  that  this  type  of  extrapolation  may  be  acceptable  where  the  depth  of  flow  over 
the  flood  plain  is  less  than  20%  of  the  channel  depth. 

When  the  depth  of  flow  over  die  flood  plain  is  a  significant  proportion  of  that  in  the  channel 
there  is  strong  inter-action  between  the  in-channel  and  ovcrbank  elements.  The  flood  plain  has 
ceased  to  be  a  reservoir  and  become  instead  the  flood  way  of  a  two-stage  channel,  quite  possibly 
carrying  the  bulk  of  the  discharge.  There  is  active  exchange  of  water  (and  sediment)  between  the 
channel  and  the  flood  plain  flows,  and  the  in-bank  element  cannot  be  considered  in  isolation.  For 
example,  recent  research  in  the  Science  and  Engineering  Research  Council's  flood  flume  facility  at 
Wallingford  has  shown  that  when  the  sense  of  rotation  of  the  main,  helical  flow  in  a  bend  may  be 
reversed  when  the  overbank  flow  depth  attains  50%  of  the  channel  depth  (Donald  Knight,  personal 
communication,  1991).  This  is  associated  with  the  meandering  channel  acting  as  a  slot  in  the  bed 
of  the  floodway  rather  than  as  a  sinuous  channel  in  its  own  right. 

Under  such  circumstances  it  is  unrealistic  to  expect  any  channel  flow  model  to  accurately 
predict  scour  pattern  or  depth  based  only  on  the  in-bank  flow  parameters.  Major  scour  will 
probably  be  associated  with  points  where  overbank  flow  plunges  back  into  the  channel,  and  this 
may  well  be  at  the  inner  bank,  or  at  some  indeterminate  point  controlled  by  the  pattern  of  flood 
plain  topography,  vegetation  or  development.  With  the  present  state  of  art  either  physical  or 
numerical  modeling  of  the  flow  in  the  entire  flood  plain/channel  system  would  be  required  to  make 
reliable  predictions  of  velocity  and  scour  patterns. 


CONCLUSIONS  AND  RECOMMENDATIONS 

This  study  has  assembled  a  data-set  for  256  bends  in  single-thread  rivers,  stream  and  flume 
channels.  The  data  have  been  examined  and  screened  and  are  believed  to  be  of  a  high  quality.  All 
the  observations  correspond  to  high  in-bank  flows  and  should  represent  the  channel  forming  flow. 
All  the  rivers  are  believed  to  be  flowing  through  alluvial  channels  in  as  much  as  the  bed  materials  are 
fully  mobile  for  the  reference  conditions.  On  this  basis  the  data  should  allow  application  of  the 
models  for  bend  flow  and  scour  prediction. 

The  results  indicate  that  considerable  errors  can  occur  in  the  predicted  scour  depths.  The 
empirical  method  produces  the  most  reliable  scour  estimates,  with  the  great  majority  of  predictions 
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being  within  +/-  25%  of  the  observed  maximum  bend  scour  depth.  Bridge's  bend  flow  model  does 
nearly  as  well  for  long  radius  bends,  but  massively  over-predicts  scour  depth  for  tight  bends.  As  a 
general  rule,  the  Bridge  method  should  only  be  used  for  bends  with  Rc/w  ratios  greater  than  4.  The 
Odgaard  model  was  the  least  satisfactory  of  the  tliree.  It  systematically  under-predicted  scour  depth 
by  about  50%  and  failed  to  account  for  the  deep  scour  holes  found  in  large  rivers.  Problems  with 
the  model  center  on  the  entrainment  function  which  does  not  predict  any  scour  in  bed  materials 
coarser  than  sand,  and  in  large  rivers  where  the  program  crashed  due  to  the  prediction  of  negative 
depths  at  the  inner  bank. 

There  was  no  evidence  that  revetted  bends  suffer  deeper  scour  than  free  meanders,  at  least  in 
terms  of  the  dimensionless  scour  depth  (dmax/^bar)  scaled  on  the  mean  depth  at  the  upstream 
crossing.  This  may  be  because  the  mean  crossing  depth  is  itself  a  function  of  bank  stiffness,  and  so 
bank  effects  are  accounted  for  in  its  use  as  a  scour  scaling  parameter.  If  true,  this  further  strengthens 
the  range  of  applicability  of  the  empirical  approach  adopted  by  Thome  (1988). 

Extrapolation  of  the  approaches  tested  here  to  higher  flows  is  limited  by  the  inter-action 
between  in-bank  and  out  of  bank  components  of  the  flow  for  overbank  events.  Where  this  is 
negligible,  notional  walls  may  be  envisaged,  separating  the  channel  flow  from  the  flood  plain  storage 
and  allowing  the  delineation  of  geometrical  and  hydraulic  parameters  for  the  in-channel  portion  of  the 
flow.  Hence,  provided  that  the  water  level  associated  with  the  design  flood  is  only  slightly  above 
the  top  of  the  revetment  and  remains  mostly  in-bank,  then  the  methods  described  here  can  be  used 
with  caution.  However,  when  the  overbank  flow  depth  is  a  significant  proportion  of  the  channel 
depth,  the  import  and  export  of  water  from  the  channel  to  the  flood  plain  is  not  negligible.  Flow 
patterns  and  bed  scour  will  be  dominated  by  water  and  momentum  fluxes  which  cannot  be  predicted 
by  these  (or  any  other)  in-bank  flow  models  and  it  would  not  be  wise  to  attempt  to  extrapolate  the 
methodology  to  such  flows.  If  significant  out-of  bank  flows  over  and  behind  the  revetment  are  to  be 
allowed  for  in  the  design  then  some  form  of  flood  plain  water  and  sediment  model  is  required. 
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Fig.  1  Dimensionless  Maximum  Bend  Scour  Depth  as  a  function  of 
Bend  Geometry  for  Data  from  Thome  and  Abt  (1990) 


Rc/w 


Fig.  2  Dimensionless  Maximum  Bend  Scour  Depth  as  a  function 
of  Bend  Geometry  for  Red  River,  1981 


Fig.  3  Dimensionless  Maximum  Bend  Scour  Depth  as  a  function 
of  Bend  Geometry  for  Red  River  1969 


Fig.  5  Dimensionless  Maximum  Bend  Scour  as  a  function  of  Bend  Geometry 
for  data  from  various  researchers  Note:  open  symbols  from  River  Ganges 
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Fig.  7  Observed  vs  Predicted  Maximum  Bend  Scour  Depths 
for  Data  from  Thome  and  Abt  1990 
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Fig.  8  Observed  vs  Predicted  Maximum  Bend  Scour  Depths  for 
Data  from  Red  River  1981  (Thome,  1988) 


•  • 

• 

• 

• 

. •_ 

• 

• 

9 

• 

1  *  •  1 

!•  ,  l*  1 

_ 1 * 

i 

1 

i 

t 

5  10  IS 


•  Bridge's  dmax 
+  Od  guard's  dmax 

•  Empirical  dmax 


Observed  Scour  (m) 


3 


Predicted  Scour  Depth  (m)  Predicted  Scour  (m) 


Fig.  9  Observed  vs  Predicted  Maximum  Bend  Scour  Depths  for 
Data  from  Red  River  1969 
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Fig.  10  Observed  vs  Predicted  Maximum  Bend  Scour  Depth 
for  British  Gravel-Bed  Rivers 
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Fig.  1 1  Observed  vs  Predicted  Maximum  Bend  Scour  Depths 
for  data  from  various  researchers 
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•  Odgaard's  dmax 

•  Empirical  dmax 
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Fig.  12  Observed  vs  Predicted  Maximum  Bend  Scour  Depths 
for  Flume  Data 
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Error  (%)  Error  (%) 


Fig.  13  Errors  in  Predicted  Maximum  Bend  Scour  Depth  for  Natural  River  data 
from  Thome  and  Abt  (1990) 


♦  Bridge  error 

♦  Odgaard  error 

♦  Empirical  error 


Fig.  14  Errors  in  Predicted  Maximum  Bend  Scour 
Depths  for  Red  River  1981  data 
(Note:  empirical  method  was  developed  from  these  data) 


• 

*. 

. 

• 

L* 

*  O  •• 

• 

• 

•• 

i 

t 

1 

y  O’ 

■  • 

• 

►  ++  +* 

•  Bridge  error 
+  Odgaard  error 

•  Empirical  error 


40 


Error  (%) 


Fig.  15  Errors  in  Predicted  Maximum  Bend  Scour  Depths  for  Red  River  1969  data 
(Note  empirical  method  was  developed  for  this  river) 
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•  Bridge  error 
+  Odgaard's  error 

•  Empirical  error 


Fig.  16  Errors  in  Predicted  Maximum  Bend  Scour  Depth 
for  Bridsh  Gravel-Bed  Rivers 
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•  Bridge  Error 

•  Empincal  error 


Fig.  17  Errors  in  Predicted  Maximum  Bend  Scour  Depths 
for  data  from  various  researchers 
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Fig.  18  Errors  in  Predicted  Scour  Depths 
for  selected  Flume  Data 
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Fig.  19  Observed  vs  Predicted  Maximum  Bend  Scour  Depth  for  All  Data 
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Observed  Maximum  Scour  (m) 


Fig.  20  Errors  in  Predicted  Maximum  Scour  Depths  for  Ail  Data 
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Rc/w 


Enor  (%)  Predicted  Scour  (m) 


Predicted  Scour  (m) 


Fig.  23  Semi-log  Plot  of  Dimensionless  Maximum  Bend  Scour  Depth 
for  Revetted  Bends 


(Rc/w)  -  2 

Fig.  24  Observed  vs  Predicted  Maximum  Scour 
Depths  for  Revetted  Bends 
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•  Revetment  dmax 


Observed  Scour  (m) 


dmax/dbar 


Fig.  25  Errors  in  Predicted  Maximum  Scour 
Depths  for  Revetted  Bends 
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Fig.  26  Semiiog  Plot  of  Dimensionless  Maximum  Bend 
Scour  Depth  for  Free  Meanders 
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TABLE  l  -  BASE  DATA  FOR  NATURAL  RIVERS  FROM  THORNE  AND  ABT  (1>W)  _ _ _ 

RESEARCHER  RIVER  SITE  BEND  RADIUS  OF  MEANDER  WIDTH  MEAN 

NUMBER  CURVATURE  WAVELENGTH  (m)  DEPTH 


r*OOfs^'0'«'0'onvnvin\Ciona 


nfscNNfvnntNNiNh 
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53'^S|cssr'0<NO<^oof^o  —  «  n  - 

n  r-irtM-.~N(NWfNsrN 


de  VrieiwyGeldorf  Dommel  The  Netherlands  2  14. JO  63.40  6.00  0.50  N 

Dietrich  &  Smith  Muddy  Creek  Wyoming  I  800  50.00  4  00  0  40  N 

A.I.  Odgaard  |e,  Nishnabotna  [low; _  1  233  00  1120  00  48  00  2  05  N 

Explanation  N  =  Natural 


TABLE  I  (Continued)  -  BASE  DATA  FOR  NATURAL  RIVERS  FROM  THORNE  AND  ABT  (1990) 
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R  =  Rough  R  =  Ripples  S  =  Straight 

I  *  Intermediate  D  =  Dunes  M  =  Meandering 

S  =  Smooth  P  =  Plane  B  =  Braided 


TABLE  1  •  CALCULATED  DATA  FOR  NATURAL  RIVERS  FROM  THORNE  AND  ABt  0990) 


RcAv 

dntx/dbir 

jTggggi 

IHsiH 

BRIDGES 

dmu 

(m) 

ODGAARD'S 

dmax 

(m) 

EMPIRICAL 

dmix 

(m) 

BRIDGES 

ERROR 

(%> 

ODGAARD'S 

ERROR 

_ 

EMPIRICAL 

ERROR 

<*> 

2.87 

208 

1.35 

1.99 

0.833 

1.36 

47.41 

-38.30 

0.98 

2.16 

1.64 

5.3522 

1.3321 

226.35 

-18.77 

0.75 

1.88 

1.73 

6.087 

1.6334 

251.85 

-547 

0.88 

1.90 

1.82 

5.6725 

2.0502 

211.68 

12.65 

1.75 

2.22 

2.89 

4.01 

2.5383 

38.75 

-12.17 

1.25 

1.79 

1.59 

3.01 

1.4438 

89.31 

•9.19 

1.27 

2.64 

1.74 

2.62 

1.2743 

50.57 

-26.76 

1.27 

2.59 

2.05 

2.7 

1.2708 

31.71 

•3801 

2.92 

1.89 

2.3 

2.56 

1.5269 

2.55 

11.30 

-3361 

10.67 

7.92 

1.37 

5.18 

5.67 

4.2234 

6.53 

9.46 

-18.47 

26  06 

6.58 

1.47 

5.27 

6.16 

4.1647 

6.39 

1689 

-2097 

21.32 

3-76 

1.37 

5.48 

9.64 

5.2184 

7.87 

75.91 

-4  77 

43.59 

1.12 

1.65 

6.34 

20.55 

4  8469 

224.13 

-23.55 

080 

1.88 

6.55 

21.3808 

6.8718 

226  42 

4.91 

7.83 

1.55 

5.27 

4.93 

3  6566 

5.92 

-645 

30  61 

12.27 

1.94 

1.63 

6 

13  15 

4.8383 

119.17 

19  36 

4.70 

1.36 

503 

7.21 

4.1945 

694 

4334 

-16  61 

38  01 

3.80 

2.31 

1.5 

285 

0.8993 

1.27 

9000 

-40  05 

-15.14 

4.84 

1.30 

1.13 

1.94 

1.2688 

1.63 

7168 

12.28 

44.13 

17.90 

1.58 

8.75 

7.12 

5.5786 

8.57 

-1863 

36.24 

-2.04 

16.23 

1.55 

8.5 

6.55 

5.532 

8.61 

-22.94 

-34.92 

1.29 

10.00 

1.33 

7.5 

7.61 

5.7044 

9.46 

1.47 

•23.94 

26  17 

9.34 

1.56 

8.25 

7.07 

5.3473 

896 

-14.30 

-35.18 

865 

8.23 

1.80 

9 

7.9 

5.0411 

8.61 

-12.22 

-43.99 

-4.31 

13.85 

1.51 

8.25 

6.99 

5.488 

8.72 

-15.27 

-3348 

5.71 

9.64 

1.42 

7.75 

7.66 

54965 

9.18 

-1.16 

-29.08 

18.40 

9.26 

1.89 

9.5 

6.72 

50652 

8.52 

-29.26 

-4668 

•10.34 

11.75 

2.14 

11.75 

7.06 

5.5147 

8.97 

39.91 

-5307 

-23.63 

21.62 

1.65 

8.5 

6.28 

5.1692 

7.75 

-26.12 

-39  19 

-885 

11.92 

1.52 

8 

6.71 

5.2897 

8.58 

-16.13 

-33.88 

7.23 

2.72 

3.30 

165 

1.04 

0.6406 

1.07 

-36.97 

•61.18 

35.39 

2.42 

3.00 

IS 

0.865 

0.6549 

1.12 

-42.33 

-56.34 

-25.46 

2.00 

2.13 

0.85 

0.69 

0.6219 

-18.82 

-26  84 

4.85 

2.10 

4.3 

4  05 

2.726 

384 

-5.81 

-3660 

-1081 
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TABLE  3  •  BASE  DATA  FOR  REP  RIVER  HYDROGRAPHIC  SURVEY.  1981  FROM  THORNE.  1988 _ _ _ 

SITE  RADIUS  OF  WIDTH  MEAN  AVERAGE  SLOPE  FRICTION  "DISCHARGE  I  MEANDER  SINUOSITY 
(River  Mile)  CURVATURE  DEPTH  VELOCITY  FACTOR  WAVELENGTH 

_ [in] _ (mj _ (m} _ (m/s) _ («  10-3) _ (cumccs) _ Jm} _ 

422.5  1219.2  304.8  4.9  1.671  0.14397  0.19837  2495  5140  1.6 
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354  1371.6  335.3  4.6  1  836  0  13258  0  12892  2832  4704  1  03 

353.5  2743.2  365.8  4  52  1.713  0  13258  0.16031  2832  2448  I W 

352.5  914.4  396.2  4  2  1.702  0.18515  02107  2832  3984  108 

351.2 _ 762 _ 350.5  4  38 _ 1845  0  18454  018641 _ 2832 _ 1680  1  49 


TABLE  3  (Conlinutd)  •  BASE  DATA  FOR  RED  RIVER  HYDROGRAPHIC  SURVEY,  1981  FROM  TIIORNE.  1988 
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TABLE  4  -  CALCULATED  DATA  FOR  RED  RIVER  11  VDROORAPiUC  SURVEY.  1981  _ _ _ 

Rc/w  dm«/dbar  MEASURED  BRIDGE'S  ODOAARD'S  EMPIRICAL  BRIDGES  ODGAARD'S  EMPIRIC  A I 
dmw  dm  a*  dm«  dmax  ERROR  ERROR  ERROR 
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4  10  1  97  9  73  IMS  4  99  9  55  14  61  48  73  187 

4  09  2  00  9  22  6  49  4  63  8  88  29  56  49  81  -3.72 

7  50  1.75  7  92  9  35  4  53  7  89  18  07  42  78  0  35 

2  31  2  18  9.14  8  36  4  23  9  63  8  55  53  72  5  40 

2  17  2  44  1068  21  81  4  42 _ 10J2 _ 104  19 _  58  57 _ _M8 


TABLE  4  (Continued)  -  CALCULATED  DATA  FOR  RED  RIVER  HYDROGRAPHIC  SURVEY,  1981 
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TABLE  5  -  BASE  DATA  FOR  RED  RIVER  IIYDROORAPII1C  SURVEY,  1969 
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TABLE  5  (Conllnutd)  -  BASE  DATA  FOR  RED  RIVER  HYDROGRAPHIC  SURVEY,  1969 
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TABLE  <  -  CALCULATED  DATA  FOR  RED  RIVER  HYDROGRAPHIC  SURVEY.  1969 
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TABLE  6  (Continued)  •  CALCULATED  DATA  FOR  RED  RIVER  HYDROGRAPHIC  SURVEY,  1969  _ _ _ 

Rc/w I  dmM/dbar  I  MEASURED  I  ODOAARDS  I  BRIDGE'S  I  EMPIRICAL  I  BRIDGE'S  I  ODGAARD'S  [  EMPIRICAL 
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TABLE  7  •  BASE  DATA  FOR  BRITISH  GRAVEL-PEP  RIVERS  FROM  »EV  AND  THORNE  (1986)  _  _  _ 

RIVER  I  RADIUS  OF  I  MEANDER  |  WIDTH  1  MEAN  IvELOCITyI  SLOPE  I  FRICTION  I  DISCHARGE!  SINUOSITY 
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TABLE  7  (Conlinutd)  •  BASE  DATA  FOR  BRITISH  GRAVEL-BED 'RIVERS  FROM  HEY  AND  THORNE  II986) _ _ _ 

RIVER  RADIUS  OF  MEANDER  |  WIDTH  MEAN  I  VELOCITY]  SLOPE  [FRICTION  |dISCHARGE|  SINUOSITY 
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TABLE  g  -CALCULATED  DATA  FOR  BRITISH  GRAVEL-BED  RIVERS  FROM  KEY  AND  THORNE  (1986) 


Rc/w 

dmax/dbar 

BRIDGE'S 

dmax 

(m) 

EMPIRICAL 

r>-  ax 

(ITl| 

BRIDGES 

ERRCR 

(%> 

EMPIRICAL 

ERRCR 

(%) 

3.55 

1.44 

2.21 

3.45 

3.04 

56.15 

37.55 

6.77 

1.78 

2.79 

5.94 

2.78 

112.79 

-0.22 

3.01 

1.76 

2.88 

8.09 

3.39 

180.82 

17.76 
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6.13 
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7.46 
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4.39 

3.50 
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21.22 

3.96 

1.61 

2  18 
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2.62 
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20.28 

4.89 

1.62 

1.89 

1.69 

2.19 

-10.43 

15.64 

1.28 
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2.08 

4.91 

136.22 

1.19 

1.27 

1.63 

6.08 

272.83 

3.01 

2.01 

2.07 

5.26 

2.13 

154.31 

2.90 

1.22 

1.53 

1.74 

5.21 

199.57 

1.42 

1.83 

1.72 

6.71 

289.90 

3.02 

1.59 

0.81 

0.89 

1.05 

10.47 

30.1 1 

3.84 

1.50 

0.99 

1.04 

1.29 

5.38 

30.28 

2.92 

1.61 

1.1 1 

1.06 

1.44 

4.81 

29.68 

4.71 

1.43 

1.16 

1.07 

1.52 

-7.95 

31  32 
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TABLE  8 (CoaliwMMy CALCULATED  DATA  FOR  BRITISH  GRAVEL-BED  RIVERS  FROM  HEY  AND  TIIORNE  (1986) 
I - rS7w  [  dmax/dbar  I  MEASURE)  I  BRIDGES  I  EMPIRICAL  I  BRIDGES  I  EMPIRICAL  | 


8. IB 
10.88 
1.21 
1.16 
1.79 
2.56 
7  85 
3.02 
3.07 

3  88 
2.B4 

3.65 
2.28 
2.47 

4  04 
3.77 
2.94 
2.62 
3.43 
3  73 
4.19 

5.66 
2.59 
2  60 
8.36 
2.11 
3.55 

4.68 
2.08 
2.62 
5.11 
6.14 
7.27 
8.81 
5.76 
2.50 
3.03 
2.50 

2.69 
0.76 
1.81 


3.14 

66.61 

5.10 

4.34 

146.87 

35.50 

3.16 

37.86 

9.65 

3.65 

31.32 

-7.60 

3.37 

144.30 

0.54 

4.91 

245.26 

55.78 

1.33 

7.36 

18.12 

1.21 

18.16 

22  96 

1.41 

102.86 

59.74 

2.01 

71.62 

223.16 

34.81 

2.40 

28.71 

2.1 1 

1.35 

104.92 

9  94 

1.68 

161  84 

13  69 

3.77 

12  38 

16  24 

4.11 

5.68 

19  09 

2  82 

•1  1.07 

10  11 

2.89 

6.87 

25.17 

4.12 

35  13 

37.74 

2.86 

86.64 

41  45 

2.37 

59.61 

28  76 

2.40 

13.81 

4  17 

3.79  1 

16.25 

28  59 

4.15 

21.47 

35.18 

3  63 

36.82 

18  74 

2.51 

1  14.44 

5  08 

2.63 

33.33 

35  48 

1.59 

16.74 

28  92 

3.83 

207.49 

40.70 

2.34 

-3.28 

7.32 

2  28 

34.36 

13  82 

2.26 

93  35 

-6.64 

1.92 

-14.19 

-20.82 

2.24 

1.41 

-13.22 

4.70 

33.01 

22.45 

4.65 

56.94 

18.82 

5.99 

41.65 

16.67 

5.62 

101.99 

28.85 

5.81 

22.92 

23.84 

2.26 

9.95 

-18.61 

2.00 

40.96 

•7.00 

1.54 

85.08 

13  94 

1.65 

1 15.12 

31  74 

0.72 

-27  40 

31.55 

1  80 

•0  83 
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TABLE  9  •  BASE  DATA  SUPPLIED  BY  OTHER  RESEARCHERS _ _ _ _ _ _ _ 

RESEARCHER  RIVER  LOCATION  BEND  RADIUS  OF  MEANDER  WIDTH  MEAN  IVELOCITYl  SLOPE  FRICTION iDISCHARGdSINUOS 

No.  CURVATURE  WAVELENGTH  DEPTH  FACTOR 
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CALCULATED  DATA  FOR  BASE  DATA  SUPPLIED  OTHER  RESEARCHERS  _ _ _ 

dmaVdbar  MEASURED  BRIDGES  ODGAARD'S  EMPIRICAL  BRIDGES  ODGAARD’S  EMPIRICAL 
dmax  dmax  dm  ax  dmax  ERROR  ERROR  ERROR 
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